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The  goals  of  the  (OPTOMA  Program  ar^  to  develop  an  Ocean 
descriptive  Predictive  System  for  4he;  study'Oand  forecasting  o£> the 
evolution  of  ocean  mesoscale  features:,  •  and  the  California  Current 
System.  The ;  attainment  of  these  goals  requires  the  establishment 

of  a  broad- base  ocean  observing  and  monitoring  system  that 
includes,  e.  g.  t>  hydrographic  research  cruises,  moored  arrays,  and 
remotely  sensed  data','A  In — -particular, '  jto  forecast  the  evolution  of 
the  oceanic  flow  field,  the  observing  system  must  include  a  means 
of  obtaining  real-time,  synoptic  maps  for  the  initialization  and 
verification  of  the  dynamichT  model (s)  used.  Because  of  their 
convenience  and  rapidity,  P3  flights  to  deploy  ^AXBTs  are  a  clear 
choice  for  frequent  mappings,  of  the  study  domains.  ' 

Since  February,  ■  1983  six  OPTOMA  missions  have  been  flown,  off- 
the  California  Coast.  Of  these,  four  flights  have  been  in  the 
OPTOMA  Northern  California  domain  and  two  flights  have  been  in 
the  Central  California  domain.  M  total  of  325  AXBTs  have  been 
successfully  deployed*  In  these  regions.  Analyses  of  these  data 
reinforce  recent  discoveries  about  the  character  of  the  California 

Current  System:  the  current  regime  is  highly  variable  in  nature 
and  4t  is  comprised  of  cool  anomalies,  mesoscale  eddies,  'squirts' 
and  jets,  current  filaments  and  fronts. 

The  AXBT  synoptic  maps  are  of  greatest  utility  to  the  OPTOMA 
Program  when  they  can  be  used  in  near  real-time.  To  fulfill  this 
need,  an  -Airborne  Digital  Data  Acquisition  System*  was  fabricated  at 
XPS.  This  system,  which  is  dbuilt  around  an  HP9816  microprocessor 
and  a  Sippican  MK9  digitizing  unit,  digitizes  the  AXBT  audio 
signal,  then  stores  the  profiles  on  diskette,  y'  The  system  was 

successfully  tested  on  the  OPTOMA13P  and  15P  flights  and  can  now 
be  considered  to  be  completely  operational. 

With  the  digital  data  acquisition  problem  solved,  the  open  issue 
is  now  how  to  exploit  fully  the  aircraft  missions.  Steps  being 
taken  in  this  direction  include  expansion  of  the  onboard  software 
library  to  include  'first-look'  data  plotting  and  objective  analysis 
routines,  expansion  to  a  suite  of  aircraft  deployable  airborne 

instruments  which  includes;  e.g.,  expendable  sound  speed  profilers 
and  expendable  CTD  profilers  (to  be  in  production  soon),  and 
extension  of  the  ADDAS  to  a  multichannel,  multi-instrument 
recording  system. 
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INTRODUCTION 


OPTOMA  Overview 

The  Ocean  Prediction  Through  Observations,  Modeling,  and  Analysis 
(OPTOMA)  Program  is  a  cooperative  research  effort  between  Harvard  University 
and  the  Naval  Postgraduate  School,  and  is  sponsored  by  the  Office  of  Naval 
Research.  The  first  phase  of  OPTOMA,  initiated  in  March  1982,  defined  the 
two  overall  goals  of  the  program:  (1)  to  develop  an  Ocean  Descriptive 

Predictive  System  (ODPS-  a  portable  set  of  dynamical  and  statistical  models, 
and  an  observational  network  designed  for  the  study  of  oceanic  mesoscale 
physics),  and  (2)  to  investigate  the  dynamics  of  the  California  Current  System 
(CCS) . 

To  attain  these  goals,  an  intensive  field  program  was  scheduled  for  the 

years  1982  through  1984.  A  total  of  fifteen  quasi- synoptic  studies  have  been 

accomplished  to  date  in  the  OPTOMA  Northern  California  (NOCAL)  and  Central 
California  (CENCAL)  study  domains.  Figure  1.  There  have  been  16  research 

cruises  and  4  research  flights  in  the  NOCAL  domain,  and  5  research  cruises 

and  2  research  flights  in  the  CENCAL  domain.  The  total  number  of  XBT's, 

CTO's,  and  AXBT's  deployed  during  Phase  I  surpasses  3000,  and  additions  to 

the  data  base  continue  on  approximately  a  bimonthly  basis. 

.  ne  prototype  ocean  prediction  experiments,  OPTOMA5  and  OPTOMA11,  were 
a  series  of  four  and  six  cruises  (and  one  overflight)  spanning  five  and  eight 
weeks  in  the  summers  of  1983  and  1984,  respectively.  These  series  were 

planned  so  as  to  use  the  quasigoostrophic  streamf  unction  evaluated  from  the 
dynamic  height  field  acquired  from  the  first  cruise  as  an  initialization  field  for 
the  Harvard/NPS  quasigeostrophic  numerical  model.  The  model  then  forecast  the 
evolution  of  the  input  field.  Subsequent  cruises  provided  verification  and  re- 


Figure  1:  The  NOCAL  and  CENCAL  subdomains  of  the  OPTOMA 
Program.  Isobaths  are  shown  in  meters. 
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initialization  fields.  In  these  first  real-time  ocean  forecasting  experiments  there 
was  encouraging  agreement  between  the  forecast  and  verification  fields 
(Robinson,  et.al,  1984).  From  these  numerical  and  statistical  analyses  of  the 
hydrographic  data,  a  new  picture  of  the  CCS  emerged,  which  includes  cool 
anomalies,  mesoscale  eddies,  ’squirts'  and  jets,  current  filaments  and  fronts 
(Mooers  and  Robinson,  1984;  Rienecker,  et.al.,  1985). 

It  is  now  recognized  among  the  OPTOMA  researchers  that  the  initial  stages 
of  the  project  have  evolved  sufficiently  to  enter  a  second  phase  (OPTOMA 

Workshop  Report,  1985) .  This  phase  will  extend  the  two  goals  of  the  first 
phase  and  add  a  third.  The  three  new  goals  are: 

(1)  the  developmental  ODPS  tested  in  Phase  I  will  be  refined  and 
implemented  as  a  prototype  system  in  the  CCS, 

(2)  the  ocean  features  revealed  in  the  'new  picture'  of  the  CCS  will  be 
examined  analytically  to  determine  the  dynamical  processes  responsible  for 
their  formation,  evolution,  and  destruction,  and 

(3)  new  applications  of  the  ODPS  to  related  areas  of  interest  (e.g  , 

biology,  chemistry,  acoustics,  optics,  marine  meteorology)  will  be  initiated. 

Role  of  AXBT  Surveys  in  the  ODPS 

Because  of  their  convenience,  speed,  near  all-weather  capability,  and  the 

greater  possibility  for  areal  expansion  of  the  study  domain,  P3  flights  to 
deploy  AXBTs  are  clearly  desireable  for  the  satisfaction  of  the  first  Phase  2 

goal:  an  operational  ODPS.  Combined  with  moored  and  drifting  buoys  equipped 
with  telemetering  current  meter  and  thermistor  sensor/chains  and  with  satellite 
remotely  sensed  data,  monthly  AXBT  flights  would  complete  an  observing  system 
which  could  provide  the  data  base  for  the  continually  functioning  ODPS.  Other 
systems;  e.g.,  acoustic  tomography,  inverted  echo  sounders  and  bottom 
pressure  gauges,  SOFAR  and  RAFOS  floats  could  play  significant  roles,  too,  or 
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even  serve  as  alternatives.  The  oceanographic  cruises  would  be  focused  on 
special  physics  and  feature  studies,  as  opposed  to  quasi- synoptic  mapping  of 
the  domain.  They  would  also  be  efficiently  integrated  with  those  cruises 
required  for  the  installation,  maintenance,  and  retrieval  of  instrumentation. 

Scope  and  Purpose  of  this  Report 

The  AXBT  flights  to  date  have  been  used  (1)  to  examine  the  effectiveness 
of  the  AXBT  as  a  data  acquisition  device  within  OPTOMA  (OPTOMA3P),  (2)  to 
obtain  a  large-scale  mapping  of  the  CCS  from  Pt.  Sur  to  Pt.  Arena 
(OPTOMA8P) ,  (3)  to  obtain  the  AXBT  data  in  a  digital  form  (OPTOMA11P) 

using  the  digital  data  acquisition  system  developed  at  Scripps  Institute  of 
Oceanography,  (4)  to  evaluate  the  degree  of  synopticity  of  shipboard 
hydrographic  data  (OPTOMA13P,  OPTOMA15P),  and  (5)  to  develop  and  test  an 
Airborne  Digital  Data  Acquisition  System  (ADDAS)  for  the  real-time  digitizing 
and  processing  of  the  AXBT  audio  signal  (OPTOMA13P,  OPTOMA15P) .  The 
purpose  of  Part  One  of  this  report  is  to  trace  the  development  and 
implementation  of  the  airborne  ocean  thermal  structure  mapping  project  through 
elaboration  of  the  first  four  points  above.  In  Part  Two,  the  design  and 
operation  of  the  ADDAS  is  described. 


PART  ONE:  THE  OPTOMA  FLIGHTS 


I.  FLIGHT  PLANNING  AND  OPERATIONAL  PRODUCTS 


Flight  Coordination 

In  response  to  a  formal  request  to  CNO  (Chief  of  Naval  Operations) 
initiated  in  October,  1982,  a  CNO  Project  number  for  P3  services  of  Navy 
active  duty  or  reserve  squadrons  on  a  not -to -interfere  basis  was  authorized  for 
the  OPTOMA  Project.  The  flights  were  scheduled  through  the  Operations 
Officers  of  the  particular  squadron  assigned  each  flight.  Subsequent  liaison 
activities,  such  as  confirmation  of  the  flight  plan,  details  of  the  mission 
schedule,  and  transportation  of  the  AXBTs  and  equipment,  were  conducted 
either  through  the  Tactical  Coordinator  of  the  assigned  crew  from  the 
operational  squadron,  or  the  Training  Officer  of  the  reserve  squadron. 
Permission  was  granted  for  civilian  observers  to  participate  on  OPTOMA11P  and 
OPTOMA  15P  after  completion  of  a  survival  swim  test,  a  low  pressure/night 
vision  certification  test,  an  emergency  procedures  training  course,  and  a 
security  clearance  investigation. 

Planning  Factors 

All  OPTOMA  flight  plans  are  flexible.  The  willingness  and  capability  to 
accept  changes  in  the  flight  schedule  (date  or  duration)  and  track  is  essential 
because  of  the  numerous  factors  which  may  alter  the  schedule  of  a  military 
aircraft.  Given  the  date  and  duration  of  the  flight,  the  planning  of  the  track 
is  straightforward  and  follows  the  logic  of  a  standard  oceanographic  cruise 
plan.  The  relevant  points  to  consider  are:  (1)  the  cruising  speeds  and 
altitudes  of  the  aircraft  en  route  and  on  station,  (2)  the  number  of  AXBTS  to 
be  deployed,  (3)  the  minimum  desired  station  spacing  between  deployments,  (4) 
the  depth  ratings  and  associated  transmission  times  of  the  AXBTs,  (5)  the 
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number  and  frequencies  of  the  AXBT  transmission  channels,  and  (6)  whether 
or  not  a  multichannel  recording  system  is  available. 

The  cruising  speed  and  aircraft  endurance  determine  the  distance  which  may 
be  covered  in  a  single  flight,  usually  about  1600  nm  in  eight  hours.  A 
typical  value  for  the  speed  en  route  is  300  knots  (air  speed) ;  while  on  station 
the  speed  can  vary  between  160  to  210  knots.  A  slower  speed  (ca.  160 
knots),  or  a  temporary  holding  pattern,  may  be  necessary  when  a  single 
channel  recording  system  or  deep  AXBTs  are  used,  so  as  to  minimize  the 
distance  covered  between  deployments.  Aircraft  altitude  is  important  because  of 
the  direct  effect  of  increased  drag  with  lower  altitude  on  fuel  consumption;  it 
also  affects  the  tolerable  minimum  distance  between  buoy  deployments  due  to 
line-of-sight  interference  on  the  same  channel.  In  general,  the  altitude 
maintained  during  an  AXBT  mission  is  lower  than  optimal  for  a  P3  aircraft; 
therefore,  the  maximum  duration  of  the  flight  may  be  fuel-limited  and  shorter 
than  a  nominal  mission  of  eight  to  ten  hours.  It  is  best  to  request  a  maximum 
capacity  fuel  load  to  allow  for  increased  fuel  consumption  due  to  lower  altitude, 
and  to  allow  for  additional  flight  time  for  replacement  of  the  AXBT  failures 
that  invariably  occur. 

Items  (1)  through  (6)  are  interrelated.  The  shallow  (305m)  and  deep 
(760m)  AXBTs  transmit  for  200  sec  and  500  sec,  respectively.  Set-up  of  the 
data  acquisition  system,  launching  of  the  buoy,  annotation  of  the  log  sheets 
and  analog  traces,  and  storage  of  the  data  takes  approximately  five  (nine) 
minutes  per  shallow  (deep)  AXBT.  Time  is  a  most  important  factor  when  only  a 
single  channel  recording  system  is  being  used.  If  a  multichannel  system  is 
available  (e.g,  a  commercial  multichannel  audio  tape  recorder  or  the  onboard  14 
channel  audio  magnetic  tape  recorder- reproducer  system),  then  alternate  channel 


AXBTs  can  be  launched  closer  together  and  the  overlapping  radio  signals 
recorded  simultaneously.  At  present,  standard  production  AXBTs  transmit  on 
one  of  three  channels:  12,  14,  and  16;  the  frequencies  of  these  channels  are 
170.5,  172.0,  173.5  MHz,  respectively.  Channels  14  and  16  can  be  used  off 

the  northern  California  coast;  channel  12  should  not  be  used  as  it  is  the 
emergency  broadcast  channel  of  the  California  Forest  Service. 

The  number  of  AXBTs  which  can  be  deployed  is  limited  by  the  capacity  of 
the  aircraft  (approximately  78)  and  the  flight  duration.  The  number  of 
successful  AXBT  profiles  is  limited  by  the  failure  rate  of  the  AXBTs  (usually 
5  to  15%,  manufacturer  and  lot  dependent)  and  occasional  operator  errors.  No 
radio  frequency  transmission,  or  a  weak  and  wavering  transmission,  are  the 
most  common  reasons  for  an  AXBT  failure.  In  rough  sea  conditions  an  AXBT 
may  be  lost,  or  short  data  dropouts  may  be  encountered,  due  to  washover  by 
an  ocean  wave . 

Operational  Products 

After  each  flight,  the  AXBT  profiles  are  edited:  failures  not  identified 
during  the  flight  are  deleted  from  the  data  set;  blunder  points  are  deleted 

from  the  profiles;  bottom  'tails'  on  the  profiles  are  truncated;  times,  dates, 
and  positions  are  corrected;  and  to  smooth  the  profiles  a  three-point  Tukey 

(median)  filter  is  applied  . 

A.  Objective  Analyses 

After  completion  of  the  data  editing  process,  the  AXBT  data  are 
interpolated  to  a  grid  using  an  inverse  distance  weighting  to  compute  a  grid 

point  value  from  the  four  closest  observed  values,  within  a  100  km  radius. 

The  gridded  values  are  then  contoured  using  a  bilinear  interpolation  routine. 


The  variables  most  commonly  displayed  in  this  manner  are  temperatures  at  a 
specified  depth,  depths  of  a  specified  isotherm,  and  dynamic  height  relative  to 
a  specified  reference  level. 

Proxy  salinity  profiles  are  estimated  from  the  AXBTs  and  the  T  vs  S 

relations  determined  from  the  CTD  data  acquired  during  the  quasi- synoptic 

cruises.  They  are  used  together  with  the  temperature  profiles  to  calculate 

dynamic  heights  relative  to  300m. 

Geostrophic  velocity  and  baroclinic  mass  transport  are  also  computed  and 
then  mapped  by  an  objective  analysis  (OA)  routine  which  uses  an  exponential 
decay  correlation  function,  with  a  2ero  crossing  of  approximately  50  km,  for 

the  computation  of  the  grid  point  values  from  the  observed  values. 

B .  Reports 

Prior  to  each  mission,  a  flight  plan,  which  includes  a  proposed  flight 

track,  is  provided  to  the  squadron  liaison  officer  for  review  and  approval.  A 
flight  report  and  a  data  report  are  produced  after  each  flight.  The  flight 

report  describes  the  itinerary  and  specifics  of  each  flight,  critiques  the 

operations,  and  provides  preliminary  analyses  of  the  data  acquired.  The  data 
report,  produced  after  final  editing  of  the  data,  describes  the  data  processing 
sequence  and  provides  maps  of  station  positions  and  numbers,  figures  of  the 
edited  AXBT  profiles,  isotherm  transects,  and  a  mean  temperature  profile.  Data 
from  the  OPTOMA3  and  OPTOMA8  flights  were  issued  as  separate  reports 
(Colton,  et.al.,  1985;  Wittmann,  et.al.,  1985);  data  from  OPTOMA11P,  13P  and 

15P  were  included  in  the  OPTOMA11,  13,  and  15  hydrographic  data  reports 
(Wittmann,  et.al.,  1985).  The  data  and  data  reports  from  all  of  the  flights 

discussed  in  this  report  have  been  archived  with  NODC. 
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II.  OPERATIONAL  AND  SCIENTIFIC  RESULTS  FROM  THE  OPTOMA  FLIGHTS 


A.  OPTOMA3 


This  mission  was  authorized  by  Commander,  Patrol  Wing  Pacific,  Wing  10 

(COMPATWNGTEN) ,  and  flown  on  10  Feb  83  by  an  active  duty  squadron 
stationed  at  Moffett  Field,  CA.  A  one-time  flight  waiver  was  issued  for  an 
NPS  student,  a  naval  officer,  to  participate  as  a  flight  observer  for  the 

purposes  of  liaison  with  the  flight  crew  and  documentation  of  a  typical  AXBT 
mission  for  use  in  planning  future  flights. 

Seventy -two  AXBT  stations  were  planned  at  a  spacing  of  28  km;  the  actual 
flight  track  and  station  positions  are  shown  in  Figure  2.  A  shallow  (305m)  or 
deep  (760m)  Sippican  AXBT  was  deployed  at  each  station.  The  aircraft 

maintained  an  altitude  of  approximately  1500  ft  and  an  airspeed  of  210  knots 
during  the  drops.  To  minimize  spatial  separation,  the  average  time  between 
stations  was  4.5  minutes,  which  resulted  in  truncation  of  some  of  the  deep 
AXBTs  which  transmit  for  7  minutes  or  more.  The  AXBT  profiles  were 
recorded  on  audio  tapes,  using  the  onboard  14-channel  audio  recorder,  and  as 
analog  traces  using  two  lofargram  recorders.  These  recorders  were  operated 
on  the  AXBT  VHF  channels  14  and  16.  Station  positions  were  obtained  from 
the  aircraft's  Inertial  Navigation  System  (LITON-72  INS)  with  hourly  updates  by 
radar  and  Tactical  Air  Navigation  (TACAN) .  For  a  description  of  the 

navigation  systems  used,  and  for  estimates  of  the  position  accuracy  from  each, 
see  Appendix  A. 

The  analog  traces  were  visually  digitized  using  a  gridded  template  overlay 
which  is  manufactured  to  represent  the  Navy  frequency-to-temperature  and  fall- 
rate  to  depth  equations.  This  procedure  yields  approximately  25  to  30  points 
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per  profile  with  accuracies  in  temperature  and  depth  of  i,0.40C  and  1,6.25m, 
respectively.  These  data  were  then  transferred  to  the  IBM  3033  mainframe 
computer  at  NPS  and  edited  by  removing  obvious  AXBT  failures  that  were  not 
identified  during  the  flight  or  digitization  procedure.  Of  the  70  AXBTs 
deployed,  eight  instrument  failures  and  one  suspected  instrument  failure  were 
edited  from  the  final  data  set,  yielding  a  retention  percentage  of  about  87%. 

The  warmest  temperatures  (15C)  in  the  OPTOMA3  sea  surface  temperature 
(SST)  field.  Figure  3,  were  in  the  southwestern  quadrant  of  the  domain,  and 
the  coolest  temperatures  (less  than  13C)  were  observed  in  the  northwestern 
quadrant  of  the  domain.  There  was  also  a  meandering  temperature  front  of 
0.8C  range,  oriented  northeast-southwest  through  the  center  of  the  domain. 
The  front  exits  from  the  western  side  of  the  domain  but  this  result  should  be 
viewed  with  caution  as  this  is  an  area  which  lacks  observations  (Fig.  3;  c.f., 
Fig.  B-l).  In  the  corresponding  surface  dynamic  height  map.  Figure  4,  there 
was  a  cyclone  of  at  least  3  dyn-cm  range  and  ca.  lOO  km  diameter  in  the 
northwestern  quadrant.  There  was  a  meandering  onshore  current,  parallel  to 


the  temperature  front 

and 

there  was 

evidence  of  southward  flow 

in 

the 

southeastern 

quadrant 

of 

the  domain 

and 

an  anticyclonic  feature 

in 

the 

northeastern 

quadrant . 

The 

overall  dynamic 

height  range  observed 

in 

this 

mapping  was  17  dyn-cm. 


B .  OPTOMA8 


After  OPTOMA3,  COMPATWNGTEN  at  Moffett  Field  received  cutbacks  in 
funded  flight  hours  and  fuel  allotments;  therefore,  dedicated  eight-hour 
OPTOMA  flights  could  not  be  scheduled.  Instead,  attempts  were  made  to  drop 


Figure  2:  Flight  track  and  station  positions  for  0PT0MA3. 
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Figure  4:  Objective  analysis  of  surface  dynamic  heigh 
relative  to  300m  (0PT0MA3) . 


a  limited  number  of  AXBTs  whenever  possible  near  the  OPTOMA  domains  during 
the  course  of  operational  missions.  Usually,  only  ten  or  so  AXBTs  could  be 
deployed  on  a  flight  and  the  positions  of  the  drops  could  not  be  predetermined 
or  guaranteed  to  be  within  the  study  domains.  Between  February  and 
December  1983,  a  total  of  72  AXBTs  were  deployed:  20  each  on  three 

consecutive  days  in  April,  south  and  west  of  the  CENCAL  domain,  and  12  on 

20  July,  north  of  the  NOCAL  domain  over  the  Mendocino  Escarpment  (Figure 

1). 

It  became  clear  that  alternative  P3  sources  would  have  to  be  sought  and  in 
November  a  request  was  placed  through  the  Commander,  Antisubmarine  Warfare 
Wing,  Pacific  (COMASWWINGPAC)  San  Diego,  regarding  the  possibility  of  flying 
with  a  reserve  patrol  wing.  The  request  coincided  in  time  with  the  planning 

of  certain  reserve  wing  exercises  scheduled  for  December  at  Moffett  Field. 

Thus,  two  flights  were  scheduled  within  five  days  of  each  other,  enabling  a 

large-scale  wintertime  mapping  of  the  CCS  to  be  obtained. 

The  flight  tracks  and  the  stations  were  executed  as  shown  in  Figure  5. 
The  flight  on  10  December,  Leg  I,  covered  the  northern  half  of  the  domains; 

the  flight  on  15  December,  Leg  II,  covered  the  southern  half  of  the  domains. 

Data  from  the  two  flights  were  combined  and  treated  as  a  single  synoptic  data 
set,  bounded  by  Pt.  Sur  in  the  south  and  Pt.  Arena  in  the  north.  This 
combination  yielded  a  large  areal  mapping  extending  roughly  420  km  alongshore 
by  350  km  crosshore.  Station  spacing  was  about  40  km  alongtrack.  During 
Leg  I  shallow  and  deep  AXBTs  were  deployed;  during  Leg  II  only  shallow 
AXBTs  were  deployed.  The  aircraft  maintained  an  altitude  of  approximately 
1500  ft  and  an  airspeed  of  230  knots.  The  data  were  recorded  onboard  on 
audio  tapes  using  the  14-channel  recorder  and  as  analog  traces  on  two 


AXBT  :  x 


lofargram  recorders,  which  monitored  the  AXBTs  on  UHF  channels  14  and  16. 
Station  positions  were  obtained  from  the  aircraft's  INS  with  hoprly  updates  by 

radar  and  TACAN. 

The  data  were  digitized,  using  a  Hewlett-Packard  digitizing  pad,  by  Naval 
Oceanographic  Office  (NAVOCEANO)  personnel  in  Bay  St.  Louis,  MS.  Only 

inflection  points  were  digitized,  giving  an  average  of  about  15  points  per 
shallow  profile  and  about  20  points  per  deep  profile.  The  data  were  provided 
for  OPTOMA  on  magnetic  tape  about  three  weeks  after  the  flight  and  were 
transferred  to  the  IBM  3033  for  editing.  From  the  Leg  I  data  set, 

approximately  97%;  i.e.,  37  deep  and  23  shallow  casts  were  recovered.  From 

the  Leg  II  data  set,  approximately  90%;  i.e.,  55  shallow  casts  were  recovered. 

From  the  SST  map.  Figure  6,  a  cool  anomaly  (less  than  14C)  extended 
offshore  in  the  central  eastern  portion  of  the  domain.  There  was  a  southward 
temperature  gradient  of  +1C  in  100  km  in  the  southern  100  km  of  the  region, 
and  a  slightly  warm  intrusion  in  the  northwest  quadrant.  In  the  surface 
dynamic  height  map.  Figure  7,  a  cyclonic  feature  was  present  in  the  region  of 

coolest  SST.  A  southward  jet  flowed  along  the  perimeter  of  the  cyclone.  An 

eastward  jet  along  the  southern  edge  of  the  domain  was  associated  with  the 
southward  temperature  gradient.  North  of  this  jet,  in  the  western  half  of  the 

domain,  was  a  large,  irregularly  shaped  anticyclone  in  the  region  of  the  warm 

intrusion.  The  overall  dynamic  height  range  in  this  map  was  20  dyn-cm. 

C.  OPTOMAUP 

Manual  digitization  of  the  OPTOMA3  and  OPTOMA8  analog  traces  was  labor 
intensive  and  produced  profiles  of  low  resolution.  In  addition,  completion  of 
the  procedure  required  several  weeks,  effectively  negating  the  real-time 


mapping  capability  of  the  AXBT  missions.  For  the  flights  to  be  of  greatest 

utility  to  OPTOMA,  it  became  necessary  to  find  a  means  of  acquiring  the  data 
digitally  in  real-time. 

It  was  learned  that  the  analog-to-digital  signal  conversion  system  developed 
by  M.  H.  Sessions  at  the  Scripps  Institution  of  Oceanography  and  used  in  the 
NORPAX  (North  Pacific  Experiment)  (Sessions,  1980)  could  be  made  available 

for  use  on  the  OPTOMA  flights.  This  system  uses  a  redundant  pair  of 

frequency  counters  to  convert  the  incoming  audio  signal  to  a  digital  frequency. 
First-order  smoothing  of  the  profile  is  accomplished  by  averaging  the  frequency 
counts  over  a  one  second  interval,  yielding  approximately  200  points  per  profile 
with  depth  and  temperature  resolution  capability  of  1.5m  and  0.03C  (  +  1  Hz), 
respectively  (with  the  Gent  equation,  to  be  explained  in  the  OPTOMA13P 
discussion) .  This  averaging  technique  does  not  impair  the  quoted  instrument 
accuracies  of  .+0.18C  in  temperature  and  2%  in  depth.  Each  point  is  then 

stored  on  a  digital  cassette  tape.  The  digital  cassette  is  copied  postflight  to 
a  magnetic  tape  and  the  data  are  then  transferred  to  a  mainframe  computer. 

For  monitoring  the  AXBT  profile,  the  system  also  produces  an  analog  trace 

and  a  hard-copy  printout  of  the  frequencies  in  real-time.  Time  and  position  are 
displayed  in  LED  readouts.  As  a  backup  measure,  a  stereo  tape  recorder 

independently  duplicates  the  incoming  audio  signal  during  the  flight  for 
postflight  playback  through  the  digitizing  portion  of  the  system  in  the 

laboratory. 

Since  the  digitizing  system  was  in  use  elsewhere  on  the  date  of  the 

OPTOMA11P  flight,  Mr.  Sessions  volunteered  to  participate  as  the  flight 
observer  and  record  the  AXBT  profiles  on  a  stereo  tape  recorder  for  post¬ 
processing  in  his  laboratory  at  Scripps.  The  processing  was  performed  within 
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Figure  8:  Flight  track  and  station  positions  for  0PT0MA11P 


ten  days  and  the  data  transferred  to  the  IBM  3033  at  NPS  for  editing. 

Sixty-eight  AXBT  stations  were  planned  for  OPTOMA11P  at  a  spacing  of  36 

km  or  52  km;  the  actual  stations  were  executed  as  shown  in  Figure  8.  An 
approximate  250  km  square  was  sampled  in  the  NOCAL  domain  with  an  along- 
track  spacing  of  about  36  km.  The  aircraft  maintained  an  altitude  of 
approximately  1500  ft  and  an  airspeed  of  about  210  knots.  Sixty- six  shallow 
AXBTs  were  deployed;  there  were  six  failures  in  the  field  and  three  profiles 

deleted  postflight  for  a  retention  percentage  of  approximately  87%.  Navigational 
difficulties  arose  on  this  flight  as  a  result  of  simultaneous  failures  of  the  INS 
and  LORAN-C.  Thus,  the  navigator  had  to  determine  positions  manually  from 
radio  fixes  received  from  shorebased  stations.  Accuracy  of  position  is 

approximated  to  be  +4  km  at  the  inshore  stations  degrading  to  ^8  km  at  the 
offshore  stations . 

From  the  SST  map,  Figure  9,  there  were  (1)  a  southwestward  oriented 
temperature  front  which  exhibited  rather  sharp  bends  in  direction,  (2)  a  cool 
feature  shoreward  of  the  front  and  centered  on  Pt  Arena,  (3)  a  second 
temperature  minimum  in  the  southwest  quadrant,  and  (4)  a  broad  area  of  fairly 
constant  warm  temperature  in  the  northwest  quadrant.  From  the  surface 
dynamic  height  map,  Figure  10,  there  was  a  predominantly  westward  jet  which 
deflected  sharply  to  the  south  in  association  with  an  anticyclone  in  the 
northwestern  quadrant.  The  southern  half  of  the  domain  was  dominated  by  a 

large  cyclonic  feature  which  contained  at  least  two  relative  minima.  There  was 
a  small  anticyclonic  feature  just  inside  the  southwestern  edge  of  the  domain. 
The  overall  dynamic  height  range  in  this  map  was  23  dyn-cm. 
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Objective  analysis  of  surface  dynamic  height 
relative  to  300m  (0PT0MA11P). 


D.  OPTOMA 13P  (27  Oct  84) 


Significant  advances  in  the  OPTOMA  AXBT  program  were  made  during  the 
course  of  planning  and  executing  flights  OPTOMA  3,  8  and  IIP.  In 

particular,  comparison  of  Figure  11(a)  with  11(b)  shows  the  typical  improvement 

in  profile  resolution  gained  from  the  use  of  the  SIO  data  digitization 
(OPTOMA11P)  vice  hand  digitization  of  the  profiles  (OPTOMA3)  (for  a  complete 
comparison,  a  profile  acquired  by  the  Airborne  Digital  Data  Acquisition  System 
(ADDAS)  is  also  shown  in  Figure  11(c)). 

With  only  35  points,  the  hand  digitized  profile,  Fig.  11(a),  is  smooth.  In 
general,  subjectivity  was  more  influential  in  determining  the  final  representation 
of  the  hand-digitized  profile;  i.e.,  the  small  inversion  at  400  m  and  the  spike 
at  the  base  of  the  mixed  layer  present  in  the  analog  trace  were  concluded  to 
be  erroneous  points,  rather  than  real  structures,  and  ignored.  Adding  an 

average  of  175  points  to  each  profile,  using  the  SIO  digitization  system, 
yielded  a  more  realistic  and  objective  representation  of  the  analog  profile  (Fig. 
11(b)).  For  example,  the  'kink'  in  the  mixed-layer  at  25  m  and  the  inversion 
at  160  m  both  were  preserved.  Smaller- scale  variability  along  the  length  of 

the  analog  profile  was  maintained  in  the  SIO  digitized  profile.  The  reason  for 
the  small  inversion  at  275  m  which  was  not  observed  in  the  analog  profile  is 
not  clear,  although  it  could  have  been  due  to  a  slightly  loose  cable  connection 
at  the  sonobuoy  interconnection  box  which  affected  some  of  the  OPTOMA  IIP 
profiles. 

The  convenience  of  data  acquisition  and  improved  resolution  available  from 
the  SIO  system  suggested  seeking  permission  for  continued  use  of  the  system 
in  the  OPTOMA13  flight.  However,  the  timely  development  of  an  AXBT 
digitizing  circuit  board  for  the  Sippican  MK9  XBT  digitizing  unit  offered  a 
second  option.  Since  the  OPTOMA  shipboard  DAS  had  been  recently  upgraded 
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to  include  additional  microprocessors  and  a  MK9  unit,  development  of  an  NPS 
airborne  DAS  became  economically  feasible.  Additional  particularly  attractive 
considerations  were  that  the  data  could  be  stoi'ed  in  digital  format  and 
processed  in  real-time,  and  that  the  data  would  be  in  the  standard  OPTOMA 

format  for  ready  use  in  existing  computer  software.  A  decision  was  made  to 

fabricate  a  system  with  the  following  attributes:  (1)  there  should  be  minimal 
capital  outlay  through  optimum  usage  of  current  project  hardware,  (2)  the 

system  should  be  self-contained,  portable,  and  have  its  own  navigation 

capability,  (3)  the  system  operation  and  software  execution  should  be 
straightforward  enough  to  permit  easy  operation  by  anyone,  (4)  there  should 
be  an  independent  backup  system  in  case  of  catastrophic  failure  of  the  ADDAS, 
and  (5)  if  possible,  the  system  should  allow  multichannel  recording.  The 
ADDAS,  which  satisfied  items  (1)  through  (4),  was  assembled  and  programmed 
in  time  for  OPTOMA13P.  Upgrades  necessary  for  multichannel  recording  (Item 

5)  are  now  being  examined  (see  Part  Two  for  details). 

The  OPTOMA  13P  mission  was  flown  by  a  reserve  squadrons  VP91  and 
SAU0617  stationed  at  Moffett  Field,  CA.  Approval  for  an  NTS  observer  to 
participate  on  the  flight  was  not  received;  therefore,  a  naval  officer,  who  had 

previous  experience  with  oceanographic  AXBT  missions,  was  recruited  to  operate 
the  ADDAS. 

The  flight  was  executed  as  planned  (see  Figure  12)  except  that  the  station 
spacing  had  to  be  increased  from  37  km  to  46  km  in  order  to  avoid 

transmission  interference  from  successively  deployed  AXBTs.  Hence,  only  55 
AXBTs  were  deployed  in  the  NOCAL  domain  instead  of  the  68  planned.  A 

deep  AXBT  was  launched  at  each  station  but  the  data  only  recorded  to  305  m. 
The  aircraft  maintained  an  altitude  of  3000  to  4000  ft  and  an  airspeed  of  170 
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Flight  track  and  station  positions  for  0PT0MA13P 


to  195  knots  during  line  '.’rops.  On  this  flight,  the  aircrait's  INS  was  used  to 
determine  station  positions. 

Of  the  55  AXBTs  deployed,  the  first  four  profiles  were  not  stored  on 
diskette  due  to  a  logical  error  in  the  printing  portion  of  the  AXBT  software. 
The  error  was  avoided  in  the  remaining  drops  by  simply  choosing  not  to  make 
a  plot  of  the  profiles.  The  four  missed  stations  were  replaced  in  the  grid  on 
the  homeward  leg  of  the  mission.  The  software  error  was  corrected  postflight 
and  it  has  not  recurred.  Seven  AXBTs  did  not  transmit  an  RF  signal  and  one 
more,  which  was  edited  postflight,  had  a  weak  wavering  signal.  The  final 
failure  was  due  to  a  pre-triggering  error  in  the  MK9  (once  the  system  is 
readied  for  launch  and  waiting  for  the  AXBT  signal,  it  is  susceptible  to 
interference  on  the  receiving  channel).  In  this  case,  local  radio  traffic  on 
channel  14  triggered  the  profiling  sequence  just  before  the  launched  AXBT 
impacted  the  water  surface.  The  profile  was  not  lost,  however,  as  it  was 
successfully  recorded  on  the  audio  cassette  tape  for  postflight  playback  through 
the  ADDAS.  In  all,  42  profiles  were  recovered,  yielding  a  retention 
percentage  of  83%. 

The  major  lesson  learned  on  this  flight  was  that  the  long  transmission  times 
and  the  associated  wide  station  spacing  of  the  deep  AXBTs  reinforced  the  need 
for  a  multichannel  digitization  system.  Reducing  the  speed  of  the  aircraft  to 
compensate  for  the  longer  times  is  not  a  satisfactory  compromise  as  it  reduces 
the  areal  coverage  of  the  survey.  However,  setting  aside  the  multichannel 
issue  (and  the  minor  software  error),  the  first  test  of  the  ADDAS  was 
successful.  Its  'easy  to  operate'  attribute  was  proven  by  the  fact  that  the 
system  was  operated  by  non-NPS  personnel  on  its  maiden  flight. 
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In  the  OA  map  of  the  sea  surface  temperature,  Figure  13,  there  was  a 
positive  offshore  gradient  in  temperature  from  an  inshore  minimum  of 
approximately  13C  to  an  offshore  maximum  of  17. 6C. 

At  least  three  pairs  of  equations  exist  for  the  conversion  of  audio 
frequency  to  temperature  and  the  calculation  of  depth  for  a  given  point  from 
the  known  fall  rate  of  the  AXBT  probe.  Bane  and  Sessions  (1984) 
intercompared  AXBT  profiles  computed  from  each  of  the  three  sets  of  equations 
(Navy,  Sippican,  Gent)  against  simultaneously  obtained  CTD  profiles.  The 
objective  was  to  determine  the  most  accurate  set  to  use  when  deploying  deep 
preproduction  AXBTs  (those  manufactured  by  Sippican  in  1980  under  contract 
number  N00163-80-C-0134) .  The  results  indicated  that  dependable  results  could 
be  obtained  by  using  the  Navy  depth  equation  and  the  Gent  temperature 
equation. 

In  OPTOMA13P,  preproduction  deep  AXBTs  were  deployed;  thus,  to  assess 
the  impact  of  the  differing  equations  on  the  resulting  temperature  fields  the 
OPTOMA13P  data  were  reprocessed  using  the  Gent  equation  and  the  Navy  depth 
equation.  At  each  station,  the  profiles  produced  from  the  two  sets  of  equations 
were  differenced  and  a  mean  difference  profile  computed.  The  Sippican 
produced  temperature  profiles  were  on  average  +0.4C  warmer  throughout  the 
water  column  than  the  profiles  produced  by  the  Gent  equation. 

On  27  October,  there  were  three  AXBT  deployments  which  were  within  2 

km  and  a  few  hours  of  XBT  casts  made  by  the  NOAA  ship  MCARTHUR.  The 
XBT  profiles  were  compared  to  the  AXBT  profiles  produced  from  the  two  sets 

of  equations.  From  this  nonideal,  small  sampling,  it  was  concluded  that  the 

Gent  equations  were  more  consistent  with  the  XBT  profiles  especially  in  the 
mixed  layer.  As  in  the  Banes  and  Sessions  study,  a  +0.4C  warm  bias  was 
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observed  in  the  Sippican  temperature  profiles.  However,  because  of  the  small 
sample  size  it  was  decided  to  continue  processing  the  AXBT  data  with  the 
Sippican  equations  until  further  test  cases  could  be  added  to  the  sample.  The 

OPTOMA15  data  set  is  now  being  examined  for  additional  XBT/AXBT  pairs.  In 
the  future,  flights  which  occur  within  a  research  cruise  period  shall  include 
plans  to  rendezvous  with  the  ship  (optimal  case)  or,  at  least,  to  deploy 
AXBTs  along  the  most  recent  cruise  track. 

The  equations  now  used  in  the  MK9  software  are: 

D  =  1.554602  *  t  -  1.634E-4  *  t 1  <  3)  ( 

T  =  -126.662  ♦  0.219954  *  F  -  1.705096E-4  *  F2  <4>, 

+  7.70534E-8  *  F1  -1.7958E-11  *  F4  +  1.73823E-15  *  F*  > 

where,  D  =  depth  in  meters,  t  =  elapsed  time  in  seconds  after  probe  release, 

T  =  temperature  in  degrees  Celsius,  and  F  =  frequency  in  Hz.  The  depth 
equation  (3)  was  adjusted  from  that  previously  used  by  Sippican,  subsequent 

to  Bane  and  Sessions  (1984).  These  adjusted  equations  were  used  to  process 
the  data  from  flights  OPTOMA  IIP,  13P,  and  15P. 

Salinity  profiles  for  the  AXBTs  were  estimated  from  the  CTD  data  acquired 
during  the  quasi- synoptic  cruise,  OPTOMA  13,  in  order  to  calculate  dynamic 
heights  relative  to  300  m.  The  dynamic  height  field.  Figure  14(a),  is 

dominated  by  a  southward  current  in  the  center  of  the  domain  and  a  large 
(100  km  diameter,  7  dyn-cm  range)  anticyclone  in  the  northwest  quadrant. 
Two  smaller  anticyclonic  features  (50  to  100  km  diameter,  2  to  4  dyn-cm)  are 
also  present:  one  near  the  center  of  the  domain  on  the  inshore  side  of  the 
current,  and  the  other  on  the  offshore  side  of  the  current  in  the  southwest 

quadrant.  This  map  is  consistent  with  the  equivalent  version  produced  from 


Distance  (km) 
(a)  OPTOMA 13P 


Distance  (km) 

(b)  OPTOMA  1 3  (22  Oct  to  2  Nov  84) 


Figure  14:  (a)  Objective  analysis  of  AXBT  surface  dynamic 

height  relative  to  300  m  (OPTOMA 1 3P) . 

(b)  Objective  analysis  of  XBT/CTD  surface  dynamic 
height  relative  to  300  m  (OPTOMA13). 

The  analysis  domains  were  adjusted  to  include  only  that  area 
common  to  both  mappings  (see  Appendix  B,  Fig.  B-4) 


shipboard  observations,  Figure  14(b).  To  compare  the  fields,  the  domain  size 
was  reduced  to  200  by  160  km  so  as  to  include  only  that  area  common  to  both 
mapppings  (see  Appendix  B,  Fig.  B-4).  All  of  the  primary  discrepancies, 
viz.  :  (1)  the  alignment  of  the  large  anticyclone,  (2)  the  absence  of  the  smaller 
anticyclone  in  the  XBT/CTD  field,  present  in  the  southwestern  quadrant  of  the 
AXBT  field,  and  (3)  the  orientation  of  flow  in  the  northeastern  quadrant  of 

the  AXBT  field  can  be  attributed  to  the  absence  of  data  from  the  shipboard 
survey  (1  and  2),  or  the  airborne  survey  (3),  when  the  contoured  fields  are 

produced  from  extrapolation  of  the  existing  data.  From  this  initial  comparison, 
it  appears  that  the  hydrographic  data  acquired  over  a  twelve  day  period  could 
be  considered  to  yield  a  fairly  valid  description  of  synoptic  fields. 

E.  OPTOMA15P  {27  Jan  85) 

This  flight  was  flown  on  27  Jan  85  by  reserve  squadrons  VP91  and 
SAU0617.  Ms  Colton  (NPS)  acted  as  the  scientific  member  of  the  crew  and 
operated  the  ADDAS . 

On  this  mission,  the  NOCAL  area  was  to  be  sampled  as  a  series  of 

concentric  squares,  shrinking  in  size  from  the  outermost  square,  140  nm  (260 
km)  on  a  side,  to  the  innermost  square,  60  nm  (111  km)  on  a  side.  Sixty- 
eight  AXBT  stations  were  planned.  The  planned  station  spacing  between 
deployments  was  not  as  regular  as  usual  due  to  the  staggered  use  of  shallow 
(305  m)  AXBTS  and  deep  (760  m)  AXBTs.  The  shallow  probes  require  at 
least  five  minutes  or  28  km  for  completion  of  data  transmission  and  data 

storage;  the  deep  probes  require  at  least  9  minutes  or  46  km.  This  plan  was 
nearly  accomplished,  Figure  15,  but  slower-than-planned  airspeed  (160  knots) 
and  the  necessity  to  return  to  base  early  to  meet  an  airlift  schedule  did  not 
permit  completion  of  the  innermost  square. 
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In  ail,  50  AXBTs  were  deployed  and  44  acceptable  profiles  recovered. 
There  were  six  failures:  five  due  to  instrument  malfunctions  resulting  in  no  RF 
transmission  and  one  due  to  a  pre-mature  triggering  of  the  MK9  (a  strong 
transmission  on  Channel  14  in  the  southeast  corner  of  the  domain  consistently 
triggered  the  MK9,  thereby  prohibiting  the  use  of  any  Channel  14  AXBTs  in 
this  region) . 

The  dominant  feature  in  the  SST  map,  Figure  16,  is  the  strong  offshore 
temperature  gradient  (approximately  1.4C  in  40  km)  observed  in  the  southwest 
corner  of  the  domain.  The  remainder  of  the  SST  field  is  essentially  flat, 
displaying  only  an  approximate  0.6C  range  (11.6  to  12. 2C)  over  roughly  90%  of 
the  domain.  The  warmer  temperatures  in  the  north  and  cooler  temperatures  in 
the  south,  near  the  front,  are  consistent  with  the  larger  scale  weak 
anticyclone  and  cyclone,  respectively,  apparent  in  the  surface  dynamic  height 
map,  Figure  17.  There  is  a  seven  dyn-cm  offshore  gradient  in  the  dynamic 
height  in  the  southwest  corner  of  the  domain.  This  gradient  is  perpendicular 
to  the  SST  front,  and  it  is  indicative  of  a  southeastward  current. 

Tables  I  and  II  summarize  the  operational  specifics  and  the  AXBT  statistics, 
respectively,  of  the  OPTOMA  flights  to  date.  Table  III  summarizes  the 
oceanographic  characteristics  observed  in,  or  computed  from,  the  AXBT 
temperature  fields. 
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TABLE  II:  OPERATIONAL  SPECIFICS  OF  THE  OPTOMA  FLIGHTS 
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TABLE  I  :  AXBT  STATISTICS 


TABLE  III:  OCEANOGRAPHIC  CHARACTERISTICS  OBSERVED  IN  OR  COMPUTED 
FROM  THE  AXBT  TEMPERATURE  FIELDS 


FLT  DATE  DOMAIN  DYN  HT  TEMP 

RANGE  RANGE 
(dyn-cm)  (deg  C) 


OPT3  10  FEB  83  CENCAL  17  12.5  to  15.5 

Major  Features:  Northeastward  (onshore)  current;  cyclone  in  NW. 


OPTS  10  DEC  83  CENCAL  20  13.0  to  15.9 

15  DEC  83  NOCAL 

Major  Features:  Inshore  weak  cyclone  with  southward  jet 
along  offshore  edge,  anticyclone  in  NW,  eastward  current 
along  bottom  edge  of  domain. 


OPTll  18  JUL  84  NOCAL  23  12.6  to  15.5 

Major  Features:  Westward  jet  deflected  to  the  south  by 
anticyclone  in  the  NW,  cyclone  in  the  SE, 


OPT  13  27  OCT  84  NOCAL  20  12.1  to  17.8 

Major  Features:  Southward  current  in  center  of  domain, 
anticyclone  in  the  NW;  two  smaller  anticyclones  inshore 
of  jet  and  one  anticyclone  offshore  in  SW  quadrant. 

OPT  15  27  JAN  85  NOCAL  10  11.6  to  14.6 


Major  Features:  1.4C  temperature  gradient  in  the  SW  corner 
of  domain,  large-scale  weak  anticyclone  in  north,  central 
domain  and  weak  cyclone  in  south,  central  domain. 


OPERATION  OF  THE  ADDAS 


III.  NPS  PROTOTYPE  ADDAS 


A .  Hardwire 

The  ADDAS  consists  of  the  following  components:  (1)  an  HP9816  MC68000 
16-bit  microprocessor  with  HP9121  3.5  inch  microfloppy  dual  disk  drives,  (2)  a 
Sippican  MK9  front-end  processor  (digitizing  unit),  (3)  an  ARNAV-R40  LORAN- 
C  unit,  (4)  a  Sony  DM6-C  tape  recorder  (monoaural,  to  be  replaced  with  a 
stereo  recorder),  (5)  a  Tektronix  Type  321A  oscilloscope  (to  be  replaced  with 
a  smaller  Tektronix-211) ,  and  an  HP  Thinkjet  printer  (to  be  replaced  with  a 
HP-7470A  two-pen  plotter).  Power  for  the  system  is  obtained  by  inverting 
28VDC  power  to  115V,  60  Hz  via  a  static  inverter  manufactured  by  Avionic 

Instruments  Inc.  (model  2A125-1A). 

B .  Software 

The  data  acquisition  software,  written  in  HP  Basic  3.0,  is  a  streamlined 
version  of  the  shipboard  XBT  data  acquisition  software.  Portions  of  the 
program  which  relate  to  the  MK9  unit  were  adapted  from  the  Sippican  MK9 
software  written  for  the  HP9845  microprocessor  (Sippican  manual,  1983) .  To 
make  the  AXBT  profiles  appear  like  the  XBT  profiles  for  input  into  the 
OPTOMA  data  processing  software,  the  audio  frequencies  are  averaged  to  make 
the  profiles  a  standard  684  points  in  length  (the  background  for  this  standard 
format  is  given  in  Appendix  C).  With  a  sampling  interval  of  0.1s,  a  shallow 
AXBT  which  transmits  for  approximately  205s  will  consist  of  about  2050  points; 
therefore,  the  averaging  interval  needed  to  reduce  this  figure  to  684  points  is 
0.3s,  (3  points).  In  the  case  of  a  deep  AXBT,  nearly  5000  frequency  values 
are  transmitted;  therefore,  an  averaging  interval  of  0.7s  (7  points)  is  used. 
The  maximum  profile  depths  obtained  for  the  shallow  AXBTs  is  305  m  and  for 


the  deep  AXBTs  it  is  760  m.  The  accuracy  of  the  temperature  measurement  is 
+0.18C  and  the  accuracy  of  the  depth  measurement  is  ±5  m,  or  2%  whichever 
is  greater. 

The  ADDAS  program  is  driven  from  a  'menu'  which  has  the  following 

options : 

1)  Plot  AXBT  -  allows  plotting  of  any  AXBT  on  any  disk 

2)  Change  time  and  date  -  (executed  once  per  flight  before  first  station) 

3)  Catalog  data  disk  -  allows  review  of  contents  of  any  disk 

4)  MK9  AXBT  -  prepares  system  for  deep  or  shallow  AXBT  deployment. 

When  an  option  is  selected,  prompts  are  displayed  which  step  the  user  through 
the  inputs  required  for  that  option.  Detailed  instructions,  and  a  program 
listing  are  provided  in  Appendix  C. 

C.  Installation  of  the  ADDAS 
l)  In  the  portable  rack 

The  ADDAS  components  are  mounted  in  a  22"  X  22"  X  30"  half-rack  as 
shown  in  Figure  18.  The  MK9  unit  is  held  in  place  by  two  rack  mounting 
brackets  and  four  flat-head  screws.  The  oscilloscope,  disk  drive,  and  LORAN 

unit  are  stacked  on  top  of  the  MK9  with  firm  packing  foam  between  them  for 

protection  and  immobility.  The  Thinkjet  printer  is  set  up  on  the  sensor  station 
counter  after  takeoff  so  that  there  are  no  obstructions  in  the  paper-feed  path. 
(These  somewhat  temporary  arrangements  were  considered  acceptable  for  the 
prototype  ADDAS  since  the  system  was  usually  dismantled  immediately  after  each 
flight  for  use  of  the  components  in  other  projects;  more  frequent  flights  will 
mandate  more  secure  installation  methods.)  The  HP9816,  keyboard  and  tape 
recorder  are  placed  on  top  of  the  rack  and  are  wedged  between  wooden  blocks 
for  lateral  stability.  Vertical  motions  of  these  components  are  restricted  by 


Figure  18:  ADAS  component  installation 
in  portable  rack. 


stretching  tubing  over  them  and  anchoring  the  tubing  to  eyebolts,  which  are 
screwed  into  the  plywood  on  top  of  the  rack. 

The  power  inverter  bottom  surface  is  coated  with  thermal  conducting  grease 
and  the  inverter  mounted  to  a  12"  X  24"  X  1/8"  aluminum  baseplate,  which 
acts  as  an  external  heat  sink.  The  assembly  is  then  bolted  to  the  rear,  right 
side  of  the  rack.  The  input  to  the  inverter  is  a  10  AWG,  15  ft  long  cable 
soldered  to  a  pin  connector  on  the  front  plate  of  the  inverter.  On  the 

airframe  end  of  the  cable  are  input  and  ground  spade-type  terminal  lugs.  The 
output  to  the  system  is  a  18  AWG,  3  ft  cable  with  a  female  plug  on  the 
system  end.  The  toggle  switch  on  the  right  side  of  the  system  control  panel 
allows  28VDC  aircraft  power  to  reach  the  inverter  only.  The  lamp  next  to 
this  switch  glows  when  the  inverter  is  operating  normally. 

The  power  strip  on  the  rear,  left  side  of  the  rack  plugs  into  the  output 

cable  of  the  inverter.  All  of  the  system  components,  except  the  tape  recorder 
which  is  battery  operated,  then  plug  into  the  strip.  For  electrical  protection, 
a  switch  on  the  strip  isolates  the  system  from  the  inverter.  Each  component 
has  its  own  power  switch,  thus  a  third  level  of  power  isolation  is  achieved. 

Six  cables  (two  per  channel)  with  test  lead  connector  pins  extend  from  the 
two  rotary  switches  labelled  MK9  and  RECORDER  (SI  and  S2  in  Figure  18)  on 
the  left  side  of  the  control  panel  to  the  sonobuoy  interconnection  box.  Any 

one  of  three  channels  may  then  be  selected  and  recorded  by  the  MK9  or 

recorder.  The  oscilloscope  receives  its  input  signal  at  the  BNC  connector  just 
to  the  right  of  these  switches  on  the  control  panel. 

Finally,  communications  between  the  HP9816,  disk  drive,  MK9  and  printer 
are  accomplished  through  IEEE488  HP  Interface  busses. 


2)  In  the  aircraft 


The  installation  and  checkout  of  the  ADDAS  in  the  aircraft  requires 
approximately  one  hour . 

a)  Mechanical  connections 

Safety  regulations  for  the  P3  require  that  all  extra  equipment  brought 

onboard  the  aircraft  must  be  securely  fastened  for  takeoff  (and  preferably,  for 
the  duration  of  the  flight).  To  satisfy  this  requirement,  the  ADDAS  is 
attached  to  the  sensor  station  seat  rails.  The  seat  rails  are  two  1  inch  high, 
1  1/8  inch  wide  I-beams  spaced  14  inches  apart,  which  run  parallel  to  the 

sensor  stations.  One-quarter  inch  holes  spaced  one  inch  apart  allow  the  seats 
mounted  on  the  tracks  to  slide  to  any  position  and  to  be  locked  into  place. 

To  mount  the  ADDAS  on  the  rails,  two  segments  of  2  X  2.5  inch  u -channel 

are  bolted  to  the  bottom  of  the  half- rack.  Fig.  18  (the  oversize  channel  was 
chosen  so  as  to  avoid  scratching  the  sliding  surfaces  of  the  rails) .  Two  1/2 
inch  slots  are  drilled  14  inches  apart  in  each  segment;  the  extra  1/4  inch 
allows  for  slight  misalignment  of  the  seat  rails.  Four  pairs  of  1/4  inch  bolts 
and  nuts,  with  2  washers  each,  are  inserted  through  the  seat  rail  holes, 
securing  the  ADDAS. 

b)  Electrical  connections 

The  28VDC  aircraft  power  supplied  to  the  ADDAS  is  that  which  is  delivered 
to  Sensor  Station  2.  To  access  this  power,  the  input  terminal  lug  is 
connected  to  terminal  A24  on  terminal  board  TB431  above  Sensor  Station  2  (see 
Figure  19  and  Plate  1);  the  ground  terminal  lug  is  connected  to  a  stud  on  the 
far  right  side  of  the  terminal  board.  The  10  amp  circuit  breaker  for  Sensor 
Station  2  is  on  the  center  electrical  circuit  breaker  panel  located  just  inside 
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and  aft  of  the  fuselage  door. 

The  six  receiver  cables  of  fhe  ADDAS  are  routed  beneath  the  triangular 

partition  aft  of  the  ADDAS,  then  up  to  the  sonobuoy  interconnection  box 

located  on  the  mid-fuselage  bulkhead.  At  present,  only  four  of  the  six  cables 
are  connected  on  each  flight;  the  remaining  two  allow  for  expansion  to  a 
third  channel  AXBT.  The  test  points  used  at  the  junction  box  are  Standard 
Audio  pins  numbered  4  through  7;  these  pins  correspond  to  channel  select 
switches  A  through  D  at  Sensor  Station  1.  To  minimize  confusion,  the 

channels  are  selected  at  the  Sensor  Station  so  that  A  through  C  on  the 
ADDAS  rotary  switches  MK9  and  RECORDER  correspond  to  A  through  C  on  the 
channel  selects;  i.e., 

Pin  4:  Plug  in  'A'  lead  from  ADDAS;  select  channel  14  at  Sensor  Station  1, 

Pin  5:  Plug  in  'B'  lead  from  ADDAS;  select  channel  16  at  Sensor  Station  1, 

channel  select  'B', 

Pin  6:  Plug  in  'C'  lead  from  ADDAS;  select  third  channel  at  Sensor  Station 

1,  channel  select  'C'. 

The  final  connection  required  is  the  antenna  cable  connection  to  the  LORAN 
unit.  The  antenna  cable  is  routed  forward  approximately  25  ft  to  the  antenna 
coupler  BNC  connector  loaded  in  the  port  electronics  bay  just  forward  of  the 
RADAR  sensor  station.  A  preamplifier  is  placed  in-line  six  inches  from  the 
coupler  to  boost  the  signal  to  the  LORAN.  For  safety  purposes,  the  cable  is 
run  overhead  of  the  stations  and  secured  wherever  possible  with  plastic  tie- 
wraps. 

E.  Operational  usage 

To  understand  better  the  operation  of  the  ADDAS,  Figure  20  illustrates  an 
AXBT  deployment.  The  following  paragraph  from  Bane  and  Sessions  (1984) 


describes  the  major  events  of  the  deployment  sequence: 


An  AXBT  package  is  generally  deployed  top-end  first  from  a  chute  in 
the  lower  rear  of  the  aircraft's  fuselage  {Plate  2).  As  the  AXBT  leaves 
the  aircraft,  the  spring-loaded  wind  flap  separates  from  the  package  body 
and  pulls  out  the  cross-type  parachute.  The  parachute  stabilizes  the 
package  in  an  upright  position  and  controls  the  speed  of  the  package 
during  descent  to  the  ocean's  surface.  After  water  impact  the  seawater 
battery  activates  and  turns  on  timing  circuitry,  which  triggers  a 
mechanism  to  inflate  the  floatation  bag.  As  this  bag  inflates  it  forces  a 
plate  out  of  the  outer  cylindrical  housing  allowing  the  parachute  and  the 
outer  housing  to  be  discarded.  A  smaller  cylinder  containing  the  battery, 
electronics,  floatation  bag,  and  XBT  probe  with  wire  and  spools  remains 
at  the  ocean's  surface.  The  VHF  transmitter  is  turned  on  and  sends  an 
unmodulated  RF  signal  to  the  aircraft  on  one  of  three  possible  VHF 
carrier  frequencies  (170.5,  172.0,  and  173.5  MHz).  Approximately  40  s 

after  battery  activation,  the  XBT  probe  is  released  from  the  small 
cylinder  and  begins  temperature  profiling.  Temperature  is  converted  to 
an  audio-range  frequency  in  the  probe;  that  frequency  is  transmitted  up 
the  hard  wire  link  to  the  surface  electronics;  then  the  RF  carrier  is 
modulated  with  the  audio  frequency.  Probe  descent  time  is  later 
converted  to  depth  and  audio  frequency  to  temperature  to  produce  a 
temperature -depth  profile  by  using  appropriate  formulae.  About  1  min 
after  completion  of  the  descent  of  the  probe,  a  current  is  sent  through 
a  resistive  heater  wire  attached  to  the  inside  of  the  floatation  bag 
causing  it  to  puncture,  thereby  allowing  the  package  to  sink. 


The  RF  modulated  signal  is  accessed  by  the  ADDAS  at  the  sonobuoy 
interconnection  box.  To  record,  digitize,  and  plot  each  AXBT  profile,  the 
ADDAS  is  operated  using  the  following  procedure: 

1)  Power  up  system  (one  time  only) 

Toggle  inverter  power  switch  on  the  right  side  of  the  control  panel  to 

'OX'.  When  the  inverter  fault  lamp  glows,  indicating  that  the  inverter  is 

functioning  normally,  turn  on  the  power  strip  switch.  Next,  turn  on  the 

individual  components  in  the  following  order:  disk  drive,  CRT,  MK9, 

oscilloscope,  LORAN,  and  after  takeoff,  th<  printer. 
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Figure  20  :  Schematic  representation  of  AXBT  deployment. 


2)  Initialize  the  HP  BASIC  operating  system  and  load  the  AXBT  program 
( one  time  only  or  after  power  failure ) 

a)  Insert  the  BASIC(SYS)  disk  in  left  disk  drive  and  the  BINARIES  disk  in 
the  right  disk  drive. 

b)  Press  the  Shift  and  Reset/Pause  keys  to  trigger  system  memory  check. 

c)  When  BASIC  Ready  appears,  remove  the  system  diskettes. 

d)  Insert  'AXBT'  diskette  in  the  left  drive. 

e)  Type  LOAD  "AXBT",  press  Execute,  then  Run. 

f)  Remove  AXBT  diskette  from  left  disk  drive. 

g)  Insert  data  diskette  into  right  drive. 

3)  Press  the  'k0'  programmable  key  to  display  the  'Menu'.  Select  option  '4' 
(MK9  AXBT)  to  prepare  the  computer  for  an  AXBT  launch  (see  Appendix  A). 
If  the  graph  and  'Clear  to  launch*  cue  do  not  appear,  cycle  the  MK9  power 
switch  to  reset  the  I/O  interfaces. 

4)  Set  rotary  switches  on  the  control  panel  to  the  appropriate  channel  positions 

for  the  next  AXBT  to  be  launched,  e.g.,  'A'  for  channel  14,  'B'  for  channel 

16. 

5)  Complete  the  station  log  sheet.  The  log  sheet  has  columns  for  entries  of 
channel  number,  station  number,  time,  latitude,  longitude,  audio  cassette 
number  and  tape  position  counter,  and  comments. 

6)  Start  the  tape  recorder. 

7)  Start  the  lofargram  recorder  (done  one  time  only). 

8)  Inform  the  crew  over  the  Intercommunications  System  (ICS)  that  all  is 
prepared  for  the  next  launch;  launch  the  AXBT. 
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9)  To  listen  for  the  start  of  the  AXBT  audio  signal,  select  the  radio  receiver 
mode  on  the  ICS  mode  switch  above  Sensor  Station  2.  About  40  sec  after  the 
beginning  of  the  signal  the  probe  should  deploy.  Watch  for  evidence  of  this 
on  the  CRT,  on  the  lofargram,  and  on  the  oscilloscope,  which  should  display  a 
sinusoidal  waveform.  Note  failures  and  weak  signals  on  the  log  sheet. 

10)  At  the  end  of  the  profile:  'PAUSE'  the  tape  recorder,  make  a  hard  copy 
of  the  profile,  store  the  data,  and  annotate  the  lofargram. 

Steps  3  through  10  are  repeated  for  each  station. 


IV.  POSSIBLE  IMPROVEMENTS  TO  THE  ADDAS 


A.  Improvements  on  the  Basis  of  Human  Factors  Analysis 

There  are  some  obvious  improvements  which  should  be  made  to  the  ADDAS 
to  make  it  more  manageable  and  portable.  A  'custom-made'  rack  or  box  needs 
to  be  designed  that  includes  shelving  or  brackets  on  which  to  mount  the 
ADDAS  components.  Carrying  handles  on  the  sides  of  the  rack  would  make 
transportation  of  the  system  easier  and  safer.  There  should  also  be  a  door  or 
covering  which  can  be  latched  onto  the  front  of  the  box  to  protect  the 
installed  components  from  accidental  damage  and  the  weather.  A  small  bracket 
on  which  to  coil  extra  cable,  mounted  on  an  inside  wall,  would  avoid  tangled 
cables.  The  four  1/4  inch  bolts  and  nuts  which  fasten  the  ADDAS  to  the  seat 
rails  should  be  replaced  with  bolts  which  have  detents  as  it  is  somewhat 
difficult  to  reach  under  the  ADDAS  to  tighten  the  nuts.  Finally,  the  LORAN-C 
panel  annunciator  light,  which  flashes  within  two  nautical  miles  of  a  waypoint, 
should  be  mounted  on  the  ADDAS  control  panel  to  warn  the  operator  of  an 
impending  station. 

B.  Expansion  to  multichannel  recording 

The  HP9816/Sippican  MK9  single  channel  system  may  be  expanded  to  a 
multichannel  system  by  two  methods.  The  first  method,  which  at  present  is 
under  development  at  Sippican,  is  to  program  the  MK9  units  with  a  selectable 
IEEE488  HP  Interface  Bus  address.  In  this  way,  two  MK9  units,  each  one 
receiving  a  single  channel,  could  be  differentiated  between  and  responded  to 
by  one  microprocessor  through  a  single  interface  connection.  The  alternate 
method  would  be  to  add  a  second  IEEE488  HP  Interface  Bus  to  the  HP9816, 
thus  allowing  a  second  MK9  (with  the  same  address)  to  be  attached  to  the 
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system.  In  both  cases,  the  HP9816  software  would  be  written  to  act  upon 

'interrupts'  from  the  active  MK9  unit.  The  second  method  could  be 
implemented  at  NPS  since  it  does  not  require  any  modifications  to  be  made  to 
the  MK9  circuit  board. 

Since  parallel  development  of  the  second  method  could  not  be  accomplished 
any  sooner  at  NPS  than  the  date  set  by  Sippican  for  completion  of  their 
upgrade  (Fall  1985),  a  stopgap  method  will  be  used  in  the  summer  1985 

OPTOMA  flight.  A  stereo  tape  recorder  will  record  all  AXBT  deployments. 
However,  to  minimize  the  station  spacing,  only  shallow  AXBTs  will  be  recorded 
and  digitized  by  the  MK9  system  in  real-time.  Deep  AXBTs  will  still  be 

deployed;  but,  since  they  transmit  for  nine  minutes,  they  will  be  recorded 
only  on  the  tape  recorder  for  digitization  by  the  ADDAS  after  the  flight.  The 
data  available  in  real-time  will  be  comparable  in  time  and  spacing  to  that 
obtained  now,  and  the  addition  of  the  deep  AXBTs  to  the  data  set  after  the 

flight  will  halve  the  station  spacing  and  improve  the  data  base. 

C.  Software 

The  software  upgrades  to  the  data  acquisition  program  which  will  be 
incoporated  in  the  near  future  are:  the  computation  and  plotting  of  sound 
speed  profiles  using  salinity  information  from  the  OPTOMA  N'OCAL  domain's  mean 
T  vs  S  curve,  and  automatic  input  of  latitude  and  longitude  to  the  11P9816 
from  the  LORAN-C  across  the  unit's  RS-232  port. 

With  the  HP9816  microprocessor  as  part  of  the  ADDAS,  numerous 
applications  programs  could  be  implemented.  These  programs  would  have  to  be 
executed  after  the  final  station  is  completed  since  the  HP9816  is  not  a  multi¬ 
tasking  machine.  However,  short  programs  that  provide  a  'first-look'  at  the 


data  could  be  run  on  the  homeward  transect;  e.g.,  an  isotherm  plotting 
routine.  To  improve  the  quality  of  the  real-time  plots  an  HP7470  two-pen 
plotter  will  be  used. 

On  the  next  OPTOMA  flight,  an  attempt  will  be  made  to  run  the  objective 
analysis  programs  to  contour  sea  surface  temperature  and  depths  of  the  8C 
isotherm  (with  an  appropriate  selection  criterion  to  define  it,  the  mixed-layer 
depth  could  also  be  contoured) .  The  production  of  each  OA  map  requires 
about  15  minutes  elapsed  time;  hence,  one  map,  possibly  two,  could  be 
produced,  provided  no  errors  are  made,  on  the  homeward  transect  from  the 
NOCAL  domain  to  Moffett  Field,  which  is  approximately  a  20  to  30  minute 
flight.  Rather  than  conflicting  with  landing  procedures,  or  the  offloading  of  the 
aircraft  after  landing,  it  is  anticipated  that  in  these  initial  attempts,  the  OA's 
will  be  computed  after  offloading  of  the  ADDAS  and  reconnection  of  the 
microprocessor  to  ground  power.  Producing  the  first-look  maps  within  a  few 
hours  of  landing  will  satisfy  the  real-time  needs  of  the  OPTOMA  program  and 
will  provide  unique  post-flight  briefing  material  for  the  squadrons.  It  is 
planned  in  the  next  few  months  to  purchase  a  compiled  BASIC,  floating-point 
accelerator  hardware/software  package  for  the  HP9816.  With  this  upgrade, 
computations  will  be  performed  an  order  of  magnitude  faster.  Thus,  an  OA 
could  conceivably  be  completed  in  about  five  minutes  (elapsed  time)  on  future 
flights,  thereby  eliminating  the  need  for  post-landing  processing. 
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V.  SUMMARY 

The  OPTOMA  Program  goals  are  to  develop  an  Ocean  Descriptive  Predictive 
System  for  the  study  and  forecasting  of  ocean  mesoscale  physics,  and  to 
investigate  the  dynamics  of  the  California  Current  System.  The  attainment  of 
these  goals  requires  the  establishment  of  a  broad-base  ocean  observing  and 
monitoring  system  that  provides  data  from  at  least  the  following:  hydrographic 
research  cruises,  moored  arrays,  and  remotely  sensed  data.  In  particular,  to 
forecast  the  evolution  of  the  oceanic  flow  field,  the  observing  system  must 
include  a  means  of  obtaining  real-time,  synoptic  maps  for  the  initialization  and 
verification  of  the  dynamical  model(s)  used.  Because  of  their  convenience  and 
rapidity,  P3  flights  to  deploy  an  AXBT  grid  are  a  clear  choice  for  frequent 
mappings  of  the  study  domains. 

Since  February,  1983  six  OPTOMA  missions  have  been  flown  off  the 
California  Coast.  Of  these,  four  flights  have  been  in  the  OPTOMA  Northern 
California  domain  and  two  flights  have  been  in  the  Central  California  domain.  A 
total  of  325  AXBTs  have  been  successfully  deployed  in  these  regions.  Analyses 
of  these  data  (Part  One)  reinforce  recent  discoveries  about  the  character  of 
the  California  Current  System:  the  current  regime  is  highly  variable  in  nature 
and  it  is  comprised  of  cool  anomalies,  mesoscale  eddies,  'squirts'  and  jets, 
current  filaments  and  fronts. 

The  AXBT  synoptic  maps  are  of  greatest  utility  to  the  OPTOMA  Program 
when  they  can  be  used  in  near  real-time.  To  fulfill  this  need,  an  Airborne 
Digital  Data  Acquisition  System  was  fabricated  at  NPS  (Part  Two) .  This 
system,  which  is  built  around  an  HP9816  microprocessor  and  a  Sippican  MK9 
digitizing  unit,  digitizes  the  AXBT  audio  signal,  then  stores  the  profiles  on 
diskette.  The  system  was  successfully  tested  on  the  OPTOMA13P  and  15P 


flights  and  can  now  be  considered  to  be  completely  operational. 

With  the  digital  data  acquisition  problem  solved,  the  open  issue  is  how  to 
exploit  fully  the  aircraft  missions.  Steps  being  taken  in  this  direction  include 
expansion  of  the  onboard  software  library  to  include  'first-look'  data  plotting 
and  objective  analysis  routines,  expansion  to  a  suite  of  aircraft -launched 
instruments  which  includes;  e.g.,  expendable  sound  speed,  current,  and  CTD 
profilers  (to  be  in  production  soon) ,  and  extension  of  the  ADDAS  to  a 
multichannel,  multi-instrument  recording  system. 
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NAVIGATION  AND  POSITION  ACCURACY 

Accurate  geographic  location  of  the  AXBT  deployments  is  imperative  for  the 
correct  interpretation  of  the  acquired  data.  To  meet  this  requirement,  (at 

least)  four  navigational  systems  are  available  on  (most  of)  the  P3  aircraft  used 
in  the  OPTOMA  Program.  These  systems  are:  1)  the  LITON-72  Inertial 

Navigation  System  (INS),  2)  Tactical  Air  Navigation  (TACAN),  3)  Long-Range 

Navigation  (LORAN),  and  4)  Radio  Detection  and  Ranging  (RADAR).  The  P3-A 
aircraft,  used  occasionally  in  the  OPTOMA  program,  does  not  have  an  INS 

unless  one  is  brought  onboard  for  a  specific  mission.  The  P3-'B  Mod',  used 
most  often  in  the  OPTOMA  program,  has  an  OMEGA  navigational  system  in 
addition  to  the  above  four  systems.  Except  for  the  inertial  navigation  system, 
all  of  the  above  systems  ai*e  radio  navigation  aids  which  determine  position  by 
receiving  and  processing  radio  wave  transmissions. 

The  TACAN  system,  developed  for  the  military,  operates  in  the  Ultra- High 
Frequency  (UHF)  band  and  obtains  continuous  indications  of  distance  and 

bearing  from  surface  transmitters  located  within  300  nm,  or  line-of-sight 
distance  from  the  aircraft  (whichever  is  shorter) .  Distance  is  computed  from 
the  round-trip  travel  time  of  radio  pulse  signals  sent  from  the  surface  station. 
The  accuracy  of  this  system  is  0(1)  km,  or  3%  of  the  distance  from  the 

beacon  (whichever  is  greater) .  In  OPTOMA's  NOCAL  domain  (centered  150  km 
from  shore)  the  error  in  range  will  be  at  least  4.5  km.  Bearing  from  the 

aircraft  to  the  station  is  provided  by  the  station  which  transmits  a  reference 
signal  and  a  variable  bearing  signal.  The  signals  are  received  at  the  aircraft, 
phase-measured,  and  converted  into  an  azimuth  bearing  in  degrees,  with  (ideal) 
0.5  deg  accuracy  (or  1.3  km  at  150  km).  An  indication  of  the  position 

accuracy  for  each  flight  can  be  obtained  by  checking  the  TACAN  ground 


position  at  a  known  location  on  the  airstrip,  before  and  after  each  mission. 

The  LORAN-C  system  operates  in  the  100  to  115  kHz  frequency  band  and 
determines  position  by  measuring  the  difference  in  time  of  reception  of 
synchronized  pulse  signals  from  two  or  more  land-based  transmitters.  Accuracy 
of  geographic  position  can  vary  between  0.2  km  to  5.0  km,  depending  upon 
the  aircraft  position  with  respect  to  the  surface  stations.  In  addition  to  the 
onboard  LORAN,  the  ADDAS  has  a  LORAN-C  unit  with  similar  accuracies,  but 
more  automatic  operation. 

On  the  OPTOMA  flights,  RADAR  is  used  at  the  inshore  stations  to  measure 
range  and  bearing  with  respect  to  known  coastal  features.  The  range  is 
determined  by  measuring  the  round-trip  travel  time  required  for  the  radio  wave 
to  reach  the  object  and  return.  The  direction  is  determined  by  the  position 
of  the  rotating  antenna  when  the  reflected  portion  of  the  wave  is  received. 
Accuracy  of  the  range  measurement  is  claimed  to  be  0(1)  km;  accuracy  of  the 
bearing  measurement  is  approximately  one-half  of  the  scanning  beamwidth. 

OMEGA  is  a  global  network  of  eight  transmitting  stations.  These  stations 
transmit  in  the  Very  Low  Frequency  (VLF)  band  so  that  the  signals  are 
receivable  to  ranges  of  thousands  of  miles.  Each  station  transmits  bursts  of 
five  different  frequencies  (in  the  10.2  to  13.6  kHz  range)  multiplexed  together 
so  that  any  frequency  is  transmitted  by  only  one  station  at  a  time.  With 
eight  stations  transmitting  at  an  approximate  1.0  sec  interval  and  a  0.2  sec 
silent  interval  between  each,  the  entire  sequence  is  repeated  every  10  sec. 
The  airborne  receiver  synchronizes  to  the  OMEGA  broadcast  pattern  and 
measures  the  phase  relationship  to  each  available  station  from  the  receiver 
location.  This  phase  measurement  is  then  provided  to  a  central  processor 
where  a  circular  distance  to  each  station  is  computed.  These  distances  are 
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combined  with  a  current  aircraft  dead  reckoning  position  to  compute  the  OMEGA 
position.  Accuracy  of  the  OMEGA  system  is  consistently  ♦/-  3.7  km, 
depending  upon  the  level  of  sophistication  of  the  aircraft  roceiver/processor  and 
variations  in  the  OMEGA  signal  due  to  geophysical  effects,  e.g,  solar  activity 

(diurnal  and  intermittent)  and  polar  cap  absorption. 

The  INS  is  the  only  navigation  aid  which  is  independent  of  land-based 
stations:  the  geographic  position  of  the  aircraft  is  determined  via  measurement 

of  the  local  gravity  vertical,  which  has  a  unique  value  at  each  point  on  the 
globe.  The  position,  attitude  and  heading  of  the  aircraft  are  computed  using 

data  from  a  two-axis  accelerometer  and  a  pair  of  two-axis  gyroscopes,  which 

provide  a  locally  level  frame  of  reference  for  the  accelerometer.  The  accuracy 
of  position,  typically  1.8  km,  degrades  over  the  duration  of  a  flight  due  to 
drift  in  the  gyroscopes.  There  is  also  a  known  84.4  minute  oscillatory 

movement  (the  Schuler  oscillation)  associated  with  the  gyroscopes.  This  periodic 
motion  contributes  a  position  error  which  increases  with  time  and  is  manifested 
as  an  increase  in  drift  rate.  The  INS  used  on  the  OPTOMA  flights  internally 

corrects  for  the  Schuler  oscillation  and  also  outputs  an  estimate  of  the  overall 

drift  rate  at  the  end  of  each  mission. 

As  a  result  of  the  time-dependence  of  the  INS  drift  error,  positions  will  be 
more  accurate  in  the  first  few  hours  of  a  flight.  If  the  drift  rate  is 

considered  to  have  a  linear  temporal  dependence,  then  each  station  position  can 
be  corrected  by  an  amount  proportional  to  the  elapsed  time  between  takeoff  and 
each  deployment.  However,  to  be  more  precise  and  to  allow  for  any 
nonlinearity  in  the  error,  it  is  best  to  take  the  following  precautions  with  the 
INS:  1)  'align'  the  INS  at  the  beginning  and  end  of  each  flight  at  a  known 
location  on  the  airstrip,  and  2)  monitor  the  drift  rate  at  regular  intervals 


during  the  flight,  and  correct  for  it  in  real-time,  or  record  the  value  for 
post-flight  correction. 

On  each  flight  the  navigators  intercompare  the  position  readings  from  all  of 
the  various  aids  available  to  arrive  at  the  best  estimate  for  each  station 
position.  These  estimates  are  based  on  1)  the  frequency  of  occurrence  of  a 
particular  reading  within  a  set  of  readings  from  all  the  navigation  aids,  2) 
past  experience  on  how  errors  in  the  instruments  occur,  and  are  manifested  in 
the  readings,  and  3)  the  operating  environment  on  a  particular  day,  e.g,  rain 
and  clouds  scatter  and  reflect  the  RADAR  transmissions,  thus  reducing  its 
accuracy  and  usefulness  as  a  navigation  aid.  It  is  likely  that  the  subjectivity 
introduced  by  an  experienced  navigator  will  improve  the  the  accuracy  of  the 
position  estimate,  not  reduce  it.  By  considering  the  readings  from  all  of  the 
instruments  and  by  performing  careful  navigation  checks  throughout  the  flight 
it  is  estimated  that  the  positions  recorded  are  within  ♦/-  3.7km  (2  nm)  of  the 
true  geographic  location. 
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Figure  B-5:  Respective  positions  of  AXBTs  to 
imposed  OA  grid  (OPTOMAl 5P) .  (Sea 
surface  temperature  scaled  by  100) 
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"AXBT"  PROGRAM  INSTRUCTIONS  AND  PROGRAM  LISTING 

The  AXBT  data  acquisition  program  "AXBT"  is  a  streamlined  version  of  the 
shipboard  data  acquisition  program,  combined  with  portions  of  the  Sippican  MK9 
software  which  were  modified  for  use  on  an  HP9816  (a  BASIC  program  listing 

follows  the  instructions  section) .  The  original  shipboard  version  of  the 
program  was  written  for  an  HP9845  microprocessor  which  had  a  four  byte  real 
variable  type.  To  conform  to  an  IEEE  standard,  the  HP9816,  used  in  the 
ADDAS,  now  has  an  eight  byte  real  variable  type.  Since  the  four  byte  (six 
digit  precision)  is  adequate  for  the  OPTOMA  project,  the  eight  byte  variable  of 
the  HP9816  is  used  only  within  the  AXBT  program;  on  output  the  variables  are 
shortened  to  four  byte  real  variables,  scaled,  and  written  to  disk  as  two  byte 
integer  variables.  This  output  format  is  consistent  with  the  shipboard  program 
and  it  eliminates  the  need  for  doubled  storage  space.  With  the  standard 
format  (described  below),  all  of  the  existing  shorebased  data  editing  and 
processing  software  can  accept  AXBT  data. 

Each  data  disk  is  formatted  with  eight  files,  each  of  which  consists  of  360, 

18  byte  records.  The  AXBT  files  contain  18  profiles  of  20  records  each;  the 

first  record  of  each  profile  contains  header  information  as  shown  in  Table  C  -1. 
At  present  the  header  information,  except  for  the  time  and  date  from  the 
internal  clock,  are  entered  manually  for  each  AXBT  deployment.  The  remaining 

19  records  of  each  profile  contain  a  total  of  684  two  byte  integer  temperature 
values.  Because  the  depth  of  each  point  is  a  computed  variable,  it  is  not 
stored  on  disk.  After  the  flight  the  data  are  transferred  from  disk  to  the 
mainframe  computer  at  NPS,  at  which  time  the  depth  profiles  are  computed  and 
stored  with  the  corresponding  temperature  in  mass  storage  data  sets. 
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ARRAY 

MEMBER 


CONTENTS 


(1)  =  Hour  ♦  Second/100 

(2)  =  Latitude  (degrees  ♦  minutes/100) 

(3)  =  Longitude  -  100  (degrees  ♦  minutes/100) 

*  (4)  =  Heading  (decimal  degrees)  -  2.5  deg 

(5)  =  Wind  Direction  (degrees,  0-540  relative  to  ships  bow) 

(6)  =  Wind  Speed  (meters/sec) 

*  (7)  =  Dewpoint  (degrees  C) 

(8)  =  Air  temperature  (degrees  C) 

*  (9)  =  2m  Sea  Surface  Temperature  (degrees  C) 

*  (10)  =  2m  Salinity  (ppt) 

*  (11)  =  Boom  Sea  Surface  Temperature  (degrees  C) 

(12)  =  Station  number 

*  (13)  =  Ship  water  velocity  (meters/sec) 

*  (14)  =  Barometric  pressure 

*  (15)  =  Dynamic  height  at  10  meters 

(16)  =  Spare 

(17)  =  Number  of  points  if  CTD 

=  XBT  type  is  XBT  (  .6=T-7,  1=T-4,  2=300m  AXBT,  4=750m  AXBT) 


(18)  =  Month  ♦  day/100  ♦  Year*100 


*  Only  on  DAS  ships 

#  Only  stored  after  hydrographic  program  is  run. 

Note:  If  radiometers  are  installed,  the  output  voltages  are  packed 
and  stored  in  15  and  16. 


Table  C  -  1:  Contents  of  eighteen  variable  AXBT  header  array. 


In  the  program  AXBT  there  are  three  subprograms  which  perform  the 
eight  byte  to  four  byte  real  variable  conversions.  The  first,  "Shinit"  (Short 
Initialize),  called  at  the  beginning  of  the  program,  generates  a  look-up  table  to 
convert  single  byte  Binary  Coded  Digits  (BCD,  the  internal  HP  format)  to 
their  equivalent  decimal  values.  The  equivalence  table  is  then  stored  as  three 
arrays:  two  for  the  eight  to  four  byte  conversions  (CONI,  CON2);  one  for  the 
four  to  eight  byte  conversions:  Comb).  These  arrays  reduce  bit  manipulations 
in  the  subprograms  "Shrt"  and  "Ushrt"  (described  below'),  since  subsequent 
conversions  can  be  performed  via  a  simple  access  by  address  of  the 

equivalence  table  array.  The  second  subprogram,  "Shrt"  (Short),  called  by 
the  data  saving  subprogram  "Save",  is  used  for  writing  the  output  files  to 
disk.  In  this  routine  the  eight  byte  variables  are  converted  into  their  four 
byte  BCD  equivalents.  Upon  return  to  "Save"  the  four  byte  real  variables  are 
scaled  and  packed  into  a  two  byte  integer  array,  then  written  to  disk. 

Finally,  the  third  subprogram  "Ushrt"  (Unshort),  called  from  the  subprogram 
"Plotstd",  is  used  to  convert  the  stored  four  byte  variables  to  their  eight  byte 
equivalents  for  use  in  the  plotting  routine  "Plot". 

The  first  actions  taken  after  loading  and  running  the  program  are  the 

assignment  of  values  to  program  constants  and  initialization  of  the  running  date 
and  time  display  (Subprogram  "Showtime").  The  running  display  will  continue 
until  the  special  function  key,  'kO',  the  "Menu"  key,  is  pressed.  The  option 
menu  is  then  displayed  from  which  the  operator  may  choose  to  (1)  plot  a 

previously  stored  AXBT,  (2)  change  the  date  and  time  on  the  internal  clock, 

(3)  catalog  r.  data  disk  to  view  its  contents,  or  (4)  prepare  the  system  to 

receive  AXBT  data.  The  following  sections  outline  the  subprograms  called  by 

each  option.  Explicit  instructions  for  the  use  of  each  option  (CRT  prompts 
and  operator  responses)  follow  each  section. 


Option  PLOT  AXBT 


Selection  of  this 

option  branches  the 

program 

to  subprogram 

"Plotstd" . 

Through  a  series  of 

prompts, 

"Plotstd" 

establishes 

which  AXBT 

is  to 

be 

plotted  from  disk. 

reads  and 

converts 

the  header  information 

from 

the 

particular  file,  and 

continues 

to  the  subprogram 

"Plot" .  In 

"Plot" 

the 

temperature -depth  graph  is  prepared  (2  to  16  deg  C  on  the  abscissa;  0  to  500 
m  (800  m  for  deep  AXBTs)  on  ordinate),  the  data  are  read  via  subprogram 
"Pxbt"  and  the  profile  is  drawn.  As  the  profile  is  drawn,  the  temperature 
and  depth  values  are  written  to  the  screen  for  operator  information.  Station 
number,  time,  date,  latitude  and  longitude  are  written  to  the  screen.  Upon 
completion  of  the  profile,  a  hard  copy  of  the  plot  from  the  screen  may  be 
generated . 


Option  '1'  Prompts  and  Responses 

CRT:  Today's  Menu 

1.  PLOT  AXBT 

2. 

3. 

4. 

Enter  Number  and  push  'Continue' 
OP:  (Enter)  1  push  'CONT' 


CRT:  Plotting  Routine  for  AXBTs 

Insert  data  disk  and  follow  prompts 
Insert  data  disk  and  press  continue 
OP:  (Insert  disk  in  right  drive),  push  CONT 


CRT:  (Displays  disk  contents)  Filenumbers  and  Filetype 

Enter  file  number  to  plot 
OP:  (Enter  a  number)  1  to  8,  push  CONT 


CRT:  (If  filenumber  greater  than  8)  Enter  filenumber 

to  plot 

OP:  (Try  again) 


CRT:  (4  (If  file  does  not  contain  AXBTs)  This  file  does 


OP: 


not  contain  AXBT  data.  Do  you  want  to  plot  another? 
(Enter  either)  Y  or  N,  push  CONT 


CRT:  (If  filenumber  and  type  are  correct) 

Enter  AXBT  entry  number 
OP:  (Enter  a  number)  1  to  18,  push  CONT 


CRT:  (After  completion  of  plot)  Hard  copy?  (Y/N) 

OP:  (Enter  either)  Y  or  N,  push  CONT 


Returns  to  'Menu,  after  plot  or  ’N’  response 


Option  '2':  Change  TIME 


Selection  of  option  '2'  branches  to  the  subprogram  "Change  Time".  This 
option  should  be  executed  at  the  beginning  of  each  flight  to  ensure  that  the 
correct  date  and  time  are  stored  on  disk.  This  subprogram  makes  no 
additional  branches  and  returns  to  the  Menu  upon  completion. 


Option  '2'  CRT  Prompts  and  Operator  Responses 


CRT: 

Today's  Menu 

1. 

2.  Change  TIME 

3. 

OP: 

4. 

Enter  Number  and  push  ' 
(Enter)  2  push  'CONT' 

Continue' 

CRT: 

Time  is  hours  and 

_  minutes 

OP: 

Enter  the  correct  time  in 
Enter  correct  time 
(For  example)  1210.56 

the 

form  HHMM.SS 

CRT: 

(If  null  entry;  i.e.,  CONT)  returns  to  'Menu 
(If  valid  entry)  Updating  the  clock 

CRT: 

OP: 

Enter  date  (YYMM.DD) 
(For  example)  8407.18 

CRT: 

The  date  and  time  have  been 

corrected 

Return  to  'Menu' 


Option  *3':  CATALOG  Data  Disk 


Selection  of  option  '3'  branches  the  program  to  subprogram  'Cat  (Catalog) . 
In  'Cat'  the  directory  of  a  given  data  disk  is  read  and  displayed  on  the 
screen  for  review. 


Option  '3'  Prompts  and  Responses 

CRT :  Today's  Menu 

1. 

2. 

3.  CATALOG  Data  Disk 

4. 

Enter  Number  and  push  'Continue' 
OP:  (Enter)  3  push  'CONT' 


CRT:  Insert  data  disk  then  press  continue 

OP:  (insert  disk  in  right  drive),  press  continue. 


CRT:  Just  a  second... 

DISK  CATALOG  displayed 

Last  AXBT  entry  _  in  File  _ 

Press  Continue  when  you're  done 
OP:  Press  CONT  to  return  to  'Menu' 


Option  '4':  Mark9  AXBT 


Selection  of  Option  '4'  branches  the  program  to  subprogram  "Mark9".  The 
first  steps  in  "Mark9"  assign  the  total  number  of  measured  points  to  be  the 
standard  684  and  set  a  special  function  key,  'k4',  as  an  'Abort'  key  which  can 
be  used  to  stop  the  profiling  sequence  and  return  to  the  'Menu'. 


The  operator  is  then  requested  to  input  the  type  of  AXBT  (shallow  or 
deep)  and  certain  header  information  in  response  to  prompts;  the  time  and  date 
from  the  internal  clock  are  automatically  stored  in  the  header  array. 
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Subprogram  "Startd"  is  branched  to  next  in  which  the  data  disk  is  prepared 
for  input. 

The  microprocessor  then  sends  a  bit  message  to  the  MI<9  unit  to  query  the 
status  of  the  AXBT  digitizing  circuit  board.  If  a  positive  response  is  received 
(the  board  is  present),  then  the  microprocessor  sends  additional  bits  to  the 
MK9  to  prepare  it  for  the  AXBT  launch.  A  negative  response  indicates  that 
the  AXBT  board  is  not  present,  a  message  is  printed  and  the  entire  program 
stops.  If  an  input/output  (I/O)  error  occurs  between  the  HP9816  and  the 
N1K9,  then  a  warning  message  is  issued,  the  I/O  interface  is  cleared,  and  the 
initial  bit  message  is  sent  from  the  HP9816  again.  The  final  preparatory  steps 
before  the  launch  set  up  the  CRT  screen  for  a  real-time  plot  of  the 
temperature-depth  profile,  via  calls  to  subprograms  "Realset"  and  "Plot". 

The  system  is  now  ready  for  the  AXBT  launch  and  will  trigger  upon 
receiving  the  audio  signal  from  the  AXBT.  The  sampling  interval  of  the  MK9 
is  0.1  seconds;  therefore,  a  shallow  AXBT  which  transmits  for  approximately 
205  seconds  will  consist  of  a  total  of  2050  points  and  a  deep  AXBT  which 
transmits  for  approximately  500  seconds  will  consist  of  5000  points.  Thus,  to 
fit  into  the  684  point  record,  the  incoming  frequencies  are  averaged  in  groups 
of  3  and  7  points  for  shallow  and  deep  AXBTs,  respectively.  The  averaged 
frequency  point  is  converted  to  temperature  and  a  depth  assigned  to  in  the 
subprogram  "At".  The  point  is  plotted  on  the  screen  and  the  temperature 
saved  in  an  array.  After  the  final  point,  the  data  are  either  stored  or 
aborted  and  a  hard-copy  plot  made  according  to  the  wishes  of  the  operator. 
The  data  file  and  directory  are  closed  and  control  is  returned  to  the  'Menu' 
subprogram. 
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Option  '4' 


CRT: 


OP: 


CRT: 

OP: 


CRT: 

OP: 


CRT: 

OP: 


CRT: 

OP: 


CRT: 

OP: 


CRT: 

OP: 


CRT: 

OP: 


Prompts  and  Responses 


Today’s  Menu 

1. 

2. 

3. 

4 .  Mark9 

Enter  Number  and  push  'Continue* 
(Enter)  4  push  'CONT' 


Enter  S  for  shallow  or  D  for  deep 
(Enter  either)  S  or  D,  push  CONT 


Enter  LATITUDE  (dd.mmm) 
(Example)  38.30 


Enter  LONGITUDE  (ddd.mmm) 
(Example)  123.30 


Enter  Station  number 
(Example)  10 


(Graph  appears)  Clear  to  Launch 
(No  response  necessary) 


Enter  S  to  store  data  or  A  to  abort  cast 
(Enter  either)  S  or  A,  push  CONT 


Enter  P  to  plot 

(Enter)  P,  push  CONT  (to  plot) 

Push  CONT  to  return  to  'Menu',  no  plot. 
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"AXBT"  PROGRAM  LISTING 
IS  !  HE -STORE  "AXBT” 


20 

30 

40 

50 

!«•*  Cata  Acoulsition  Program  for  AXBT  »»» 

!•»•  to  run  on  HP21G  with  basic  3.0  *•' 

!•**  Paul  Wittmann  liav  3  .  1005  «** 

GO 

70 

! 

00 

OPTION  BASE  1 

00 

0IM  0ttt40I.CoJC3I.In 

buf fer*C2I 

100 

INTEGER  Out ( 10.2 ) 

110 

DIM  Ts<  2010  )  .Dpt  <  2010 ) 

!  temperature  and  deoth 

120 

DIM  Dir( 200 ) . Ad( 10) 

i  directory  and  header 

130 

INTEGER  I . J ,L. Iad( 36  ) , 

Datal .Datah.Data 

140  Irtit:  !•**  Initialize  variables  ««»•»»*»••  •••*#•• 

150 

Bx=l . 524 

!  drop  coeffiecant 

1G0 

Nds=GC4 

!  number  of  points 

170 

Degt*CHR3< 170) 

t  Decirees  symbol 

100 

L l$=CHRt( 10>&CHRt< 13) 

(  line  feed 

130 

L2t=Ll$&Llt 

!  double  space 

200 

Pa$=CHR$( 12  > 

!  new  oaoe 

210 

Cot”"7??" 

!  data  type 

220 

Yr*85 

!  year 

230 

Date  =  Yr»lC0t-5.2G 

•  date 

240 

Ho  i  b=7 

!  interface  adress 

250 

Xbt=720 

!  Ml<9  adress 

260 

Msut*“ ; HPC230X ,704 . 1 " 

!  3  1/2”  disk  drive 

270 

!  Msut*" = INTERNAL" 

!  5  1/4"  internal  disk  dr 

200 

riot=70G 

!  plotter  address 

290 

Prth*714 

!  printer  address 

300 

PRINTER  IS  CRT 

310 

GCSUB  Shinit 

•Set  up  Lookup  Table 

330 

Mo=INT(0ate  3500 ) 

340 

Da v=CR0UND<  FRACT<  Date 

•  100,2  ) 

350 

GCLCAR 

360 

PRINT  CHRt(lZ) 

370 

DISP  "Press  Menu  KEY  to  aet  the  menu..." 

300 

ON  KEY  5  LABEL  "Menu" 

GOTO  Menu 

330 

ON  KEY  0  LABEL  "Menu" 

GOTO  Menu 

400 

Som-TIMCDATC  MOD  8G40G 

!  sec  oast  midnight 

410 

Mour=INT ( Sgm/3G00 ) 

420 

Min«INT< (Son  MOD  3G00)/G0> 

430 

See* IN T ( Son -3600 ‘Hour 

60»Mi n) 

440 

PRINT  TA0XY( 3.15); 

450 

PRINT  "Greenwich  Mean 

Time 

460 

PRINT  DR0UN0(Hour.2>; " 

: " ;DR0UND<Mln,2 >; 

470 

PRINT  "  " ; 0R0UN0( Sec , 

2  ) 

480 

PRINT  * 

430  PRINT  *  Todays  date  — >  " ; Mo; “ / " ; Day ; * / " ; Yr 

S00  WAIT  1 

510  GOTO  330  !  Continue 

520  . . «•*•«»• 

S30  Menu:  !*•******•**••* 

540  |  •*»•»•••»•» . •  •••* 

5S0  OFF  KEY  0 

SG0  OFF  KEY  S 

570  GCLEAR 

580  PRINT  Pa* 

590  PRINT  LI*."  Today’s  MENU:'.L2* 

G00  PRINT  TAB<  5  > . " 1 .  PLOT  AXBT  " 

610  PRINT  TAB< 5  > .  “2 .  Chanqe  TIME" 

G20  PRINT  TAB(5 ), "3.  CATALOG  Data  Disk  “ 

G30  PRINT  TAB(5  ), "4.  Mark9  AXBT" 

640  PRINT  L2* 

650  DI5P  'Enter  Number  &  Push  ’CONTINUE’...  ” ; 

660  INPUT  "  '.Choice* 

670  PRINT  "Option  *";Choice* 

G80  SELECT  Choice* 

690  CASE  "1”  !  plottina  routine 

700  Rtp$="»»»"  !  real  time  plot  flap 

710  Justc»0  1  flaq  for  data  iust  collected 

720  GOSUB  Plotstd 

730  ASSIGN  @Fat  TO  » 

740  GOTO  Show  time 

750  CASE  "2"  t  reset  the  clock 

760  GOSUB  Chanae  time 

770  GOTO  Show  time 

7C0  CASE  "3-  I  cataloq  disk 

790  GOSUB  Cat 

800  GOTO  Show  time 

010  CASE  "4”  !  axbt  data  acquisition  loop 

820  Co$="XBT“ 

030  GOSUB  Mark9 

840  CASE  ELSE 

050  GOTO  Show  time 

060  END  SELECT 

070  . . ••••»•«••« . . 

880  Plotstd  . . . . . 

030  . . * . . . 

900  PRINT  Pa*. 'TlottinQ  routine  for  AXBTs" 

910  PRINT  "Insert  data  disk  and  follow  prompts" 

920  INPUT  "Insert  disk  and  oress  cont  ir.ue" ,  Y* 

330  GOSUB  Cat 

940  INPUT  "Enter  file  number  to  plot  ",X 

950  Jb*X«20 

360  IF  X>3  THEN  940 

970  ASSIGN  8Dir  TO  "DIR"SMsu* 

980  ENTER  «Dir;0ir( »  ) 

990  Y*0ir <  Jb+5  ) 

1000  IF  (Y;0>  AND  ( YC4  >  THEN  10G0 
1010  PRINT  "This  file  does  not  have"; 

1020  PRINT  " AX0T  data" 

1030  INPUT  "Oo  you  want  to  plot  another?" . Y$ 


1040  IF  YSO'Y"  THEN  RETURN 

1050  GOTO  940 

1060  ASSIGN  8Fat  TO  "FAT"&UAL* < X >&MsuS 

1070  INPUT  “Enter  AX8T  entry  number  ",En 

1000  En*En-l 

1090  ENTER  @Fat , 1+En»20 ; Out ( •  )  !  read  header 

1100  GOSUD  Ushrt 

1130  . . * . . . . 

1140  Plot:  !•■•«>••■  PLOTTING . •••*••••••••« 

1150  . * . •  •••* . 

1160  PRINT  Pa$ 

1170  Omin=0  !  Default  depths 

1180  Dmax=500 

1100  IF  Ad( 17  )=4 'THEN  OmaA=800  I  for  deeo  axbts 

1200  Tmln*2  1  temp  minimum 

1210  50*16  !  temp  span 

1220  Crt$="Y“  !  print  on  crt 

1230  GINIT  !  initialize  ut  aultiuo 

1240  GRAPHICS  ON 

1250  TmaA=Tmin+Sp 
1260  T  f ac  t  =Sp/ 100 

1270  Of ac t  =  ( Dmax -Dmin  1/ 100 

1280  DEG 

1200  VIEWPORT  0,70.5.95 

1300  CSIZE  2.5 

1310  Xmn*Tmin~8»Tfact  !  scalincj 

1320  XmA=Tmax +2*Tfact  !  factors 

1330  Ymx=0min-5»Dfact  ■  for 

1340  Ymn=Dmaxt 10*Df act  I  plot 

1350  WINDOW  Xmn.XmA.Ymn.YmA 
1360  Ntic-8 

1370  PEN  1 

1300  CLIP  Tmin.TmaA.DmaA.Dmin 

1300  AXES  Sp/2/Nt ic . < Oma* -Omin )/ 10 . Tmin.Dmin, 2 . 2 

1400  clip  orr 

1410  G0SUB  Dlab  !  Label  Depth  axis 

142C  GOSUD  Tlab  i  Label  Temperature  axis 

1430  IF  RtpI="T"  THEN  RETURN 
1440  PEN  0 

1450  G0SU0  Pxbt  !  read  axbt  data 

1460  PEN  1 

1400  FOR  1*1  TO  Nos 

1500  IF  OpUIKDmin  OR  Opt <  I  >>0ma*  THEN  1520 

1510  PLOT  Ts( I  ) , Do t ( I  ),1 

1520  NEXT  I 

1530  PENUr 

1530  IF  Ad(13)il00  THEN  Ad< 18 )*Ad< 10 >+Yr» 100 

1539  G»r=IMT( Ad( 181/100) 

1540  0at*Ad( 18  )-100*Gvr 

1541  Mo*INT ( Dat  ) 

1542  Da  *  *  INT(  1 00 » TRACT  ( Dat  )  >H 

1543  Year»INT(Ad(10 1/100) 

1555  PRINT  TABXY( 15.25 ); "Station; " ;Ad( 12  "  GMT : * ; Ad( 1  ) ; *  Date:  * ; lio ; ”/“ ;Dav 

1557  PRINT  TAG( 15  ) ; “  Lat : " ; Ad< 2 > ; "  Lono = " ; Ad< 3  )+ 100 
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1558  Y$*"  " 

1559  PRINT  TAB( 15  ) : “Hard  Coov’TY/N)* 

1560  INPUT  Y* 

15S1  IF  Y$<> “Y“  THEN  1570 

1562  ON  TIMEOUT  Hplb.l  GOTO  1558 

1563  DUMP  DEVICE  IS  Pr th 

1564  DUMP  GRAPHICS 

1565  PRINTER  IS  Prth 

1566  PRINT  "Station: ” ;Ad< 12 >; *  GMT : " ; Ad( 1 > ; "  Date: " ;Mo; ;Dav; i Yr 

1567  PRINT  ”  Lat:";Ad<2>;“  Loncj: " ; Ad( 3 H 100 

1568  PRINTER  IS  1 

1569  OFF  TIMCOUT  Hpib 

1570  GCLEAR 

1571  ASSIGN  @Fat  TO  * 

1572  RETURN 

1573  . . . . . »«*•»•*«. . . 

1574  Pxb t :  !*.«.  read  data  from  disk  .................... 

1575  !  . . . . . . . . . 

1576  L=0  !  counter 

1577  Ba*1.5240  I  drop  coeffiecent 

1570  Tstep*.3  I  tine  increment  for  shallow 

1579  IF  Ad(17)*4  THEN  Tsteo*.7  !  time  increment  for  deep 

1500  FOR  1*2  TO  20 

1501  ENTER  ®Fat,I»En»20; Iadf » > 

1502  TOR  J*1  TO  3G 

1533  L=L+1 

1534  T ime=Tstep*L 

1505  Ts( L  )  =  Iad( J  )/ ICO 

1586  Dpt ( L  )=Bx»Time  !  depth  ecjuation 

1587  NEXT  J 

1580  NEXT  I 

1589  RETURN 

1590  I . * . . 

1591  Oiab:  !••»»•  Label  vertical  "Z“-axis  ....... 

1592  . . . . . . 

1593  LOIR  90 

1594  L0RG  4 

1595  FOR  Dd*0min  TO  Dmax  STEP  < Dmax-Dmi n  )/5 

1596  MOVE  Tmin  -  .5*T fact ,Dd 

1597  LABEL  USING  "K" ;0d 
1590  NEXT  Dd 

1599  MOVE  Tmin-3.5* T fact ,Dmin  +  (Dmax  Dmin  )/2 

1600  OtS="DEPTH  (meters  )* 

1601  LABEL  USING  "K";0t* 

1602  PENUP 

1603  RETURN 

1604  I . . . . 

1605  Tlab:  I*.*  Label  temperature  axis  ......... 

1606  . . * . . 

1607  LOIR  0 

1608  LORG  6 

1609  FOR  Tp-Tmin  TO  Tmax  STEP  < Tmax-Tmin  )/Ntic 

1610  HOVE  To ,  Dmn-Ofac  t  *2 . 5 

1611  LABEL  USING  "K“;Tp 


I 
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1G 12  NEXT  To 

1813  Ot*-" TEMPERATURE  (Oaa  C>" 

1814  MOVE  Tmin+<  Tmex-Tmin  >/2 .Omin-4 .5*Ofact 

1615  LABEL  USING  “K";Ot* 

1616  PENUP 

1617  RETURN 

1618  !  . . . . * . . . 

1619  Chanoe  time: ! ••••••••••••*•••••••••••••••••* 

1620  . . . . . . . . . 

1621  PRINT  Pa$ 

1622  PRINT  "Tina  is  “;Hour;“  hours  and  “;Min;"  minutes  " 

1623  PRINT  L2* 

1624  PRINT  "Enter  the  correct  time  in  the  form  :  HHMM.SS  * 

1625  Time—999 

1626  INPUT  "Enter  correct  time  “.Time 

1627  IF  T ime*-999  THEN  RETURN 

1628  DISP  “Uodatina  the  clock  * 

1629  II* INT <  Time/ 100  ) 

1630  M* INT ( T ime - 100  *H ) 

1631  5*<Time  INK  Time  )  >'100 

1632  New*3G00'Hf60*M+S 

1633  SET  TIME  New 

1G34  INPUT  "Enter  date  ( YYMM.DO  )" .Date 

1635  UAIT  1 

1G3G  PRINT  LIS. "The  date  and  tine  have  been  corrected’ 

1637  WAIT  1 
1633  RETURN 

1033  ■  • . . . 

1640  Cat:  laaaaa  Disk  CatalOQ  ••  I  • • a  a aa a  a  a  a  a 

1641  !  »  . . . . . . a . 

1642  1  N0TC:  This  routine  loads  Oir(  I  ) 

1643  IT  Choices w "7"  THEN  GOTO  16S4 

1644  PRINT  CHRS ( 1 2  ' , T ABX Y ( 5 . 6  ) ; 

1645  PRINT  “ . aa . a . . 

1646  PRINT  TAOXYI 13.8 ), "DISK  CATALOG  OPTION- 

1647  PRINT  T ABX Y (5,10): 

1648  PRINT  . . . . . . 

1649  PRINT  TABXY(0.15)i "Follow  Dromots  at  the 

1650  PRINT  "bottom  of  the  screen.” 

1651  DISP  "INSERT  Data  Oisk,  then  oress 

1652  OISP  "•CONTINUE’." 

1653  PAUSE 

1654  DISP  "Just  a  second..." 

1655  PRINT  CHR$< 12  ),L2S, "  OISK  CATALOGUE  =’ 

1656  ASSIGN  @Dir  TO  "DIR"&MsuS 

1657  ENTER  @Dir ;Dir< •  ) 

1658  PRINT 

1659  FOR  1-1  TO  8 

1660  X-Dlr(  I»20«5M  1 

1661  rs  =  UAL$(Dir(20»I  )U .E-10) 

1662  PRINT  T  AB<  5  ) . "F ILE  “;I;H  "SDt*<X>. 

1663  IF  Dir(It2O)*0  THEN  1CG5 

1664  PRINT  T$l4,5I8i"/’&rC(6,7]&n/*&TS[l(2]&*  "GVAL*(Dir(20»I  +  l )  > 

1665  PRINT 
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1666 

1667 

1668 

1669 

1670 

1671 

1672 

1673 

1674 

1675 

1676 

1677 

1678 

1679 

1680 
1681 
1682 

1683 

1684 

1685 

1686 

1687 

1688 

1689 

1690 

1691 

1692 

1693 

1694 

1695 

1696 

1697 

1698 

1699 

1700 

1701 

1702 

1703 

1704 

1705 

1706 

1707 
1700 

1709 

1710 

1711 

1712 

1713 

1714 

1715 

1716 

1717 

1718 
1713 
1720  ! 


NEXT  I 
Xf-Dir<  3 ) 

Xe-0ir<4  ) 

IF  Xf-0  THEN  GOTO  1671 
PRINT  L2$ , “LAST  XBT  ENTRY 
ASSIGN  0Oir  TO  » 

IF  Cholce$0“3"  THEN  RETURN 

OISP  "Press  'CONTINUE*  when  you're  done.* 

PAUSE 

RETURN 


';Xes*  IN  FILE  " s Xf 


i 


Mark9-.  data  acquisition 


GOTO  1760 
GOTO  1760 


! . * . . 

Npts-684 

ON  KEY  4  LABEL  "ABORT' 

ON  KEY  9  LABEL  "CAST" 

Id*-"S" 

ASSIGN  0Xbt  TO  Xbt; FORMAT  OFF 
INPUT  "Enter  S  for  shallow  or  0  for  deep‘,Id$ 

INFUT  "ENTER  LATITUOE  ( dd.mmm  )" .Ad< 2  > 

INPUT  "ENTER  LONGITUDE  (  ddd.mmm >" .Long 
Ad(  3 )-Long- 100 

INPUT  “ENTER  STATION  NUMBER  ” ,Ad< 12 > 

Spm-TIMEOATE  MOO  86400 
Hr-INT<Spm/3G00) 

Min*INT( <Spm  MOD  3600  )/60) 

Ad<  1  )-Hr» 100*Mi n  (TIME  HIIMM 
Ad< 10  )=Uav  ‘DATE  YYMM.00 

Ad<17)«2  f  identification  for 

IT  Idt=“D"  THEN  Ad<17)»4  !  identification  for 

G05UB  Star t d  l  initialize  disk 

!  start  I/O  with  MKO  Are  the  correct  boards  present? 

OUTPUT  Xbt  USING  "*,0";05 
ENTER  Xbi  USING  "*,0" ;Resp. Board  1, Board  2 
PRINT  "BOARDS" .Resp.Ooard  1. Board  2 
IT  Resp-42  THCN  1707 
BEEP 

PRINT  "WARNING!!!  OAO  EDIU  RESPONSC* 

TRINT  "SHALL  ATTEMPT  TO  RESTART" 

WAIT  2 

GOTO  1698  !  OUTPUT  05 

IF  Board  1*3  OR  Board  2-3  THEN  1711  !  boards  are  present 

PRINT  "WARNING'!!  AXBT  BOARD  NOT  PRCSCNT" 

BECP 


Shallow 

deep 


OUTPUT  Xbt 
IT  Board  l1 
IT  Board  2' 
OUTPUT  Xbt 
OUTPUT  Xbt 
DI5P  TA0< 15) 


USING  "f.Q";l  !  boards  OK 

3  THEN  Board  mask-49 

0  THCN  Board  mask-50 

USING  "t,D“ ; Board  mask .51 ,21 . 1,25.15 

USING  "*,O“;0 

Clear  to  Launch" 

USING  "i,0" ; 1 


OUTPUT  Xbt 
PRINT  Pa® 

GOSUO  Realset  •  set  up  real  time  plot 
begin  data  acquisition  loop 
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Int*3  !  everaaina  interval  for  shallow 

IF  Id$-"0"  THEN  lnt-7  I  averaoina  intarval  for  deep 
DEEP  I  start  tfklno  data 

FOR  J-l  TO  Nots 
Data-0 

FOR  1-1  TO  Int 

ENTER  Xbt  USING  *#,2A";In  buffer* 

Oatal-NUMUn  bufferJCl]) 

Datah-NUM< In  buffer*I2I> 

Data*Datah*25G<0ataH0ata 
Time-Timet. 1  I  increment  tine 

NCXT  I 

Data*Data/ Int  !  averaoa  over  int 

GOSUO  At  !  convert  frea  to  temo 

ORAU  Teno.Oeoth  I  real  tine  olot 

Ts(J)“Temo  t  assian  to  tenD  array 

OISP  TAB( 10  >; "Tamo  "; DR0UND< Temp , 4 >, "Depth  " iDR0UND(Depth,4  ) 
NEXT  J 

DISP  "Finish  takina  data" 

INPUT  "Enter  S  to  store  data  or  A  to  abort  cast*,YS 
IF  Y**“ A"  THEN  17G0 
GOSUO  Save 

ON  TIMEOUT  Hoib, 1  GOTO  1754 
INPUT  "ENTER  P  FOR  PLOT " , Y* 

IF  YJO’P"  THEN  1752 
DUMP  OEUICE  IS  Prth 
DUMP  GRAPHICS 
PRINTER  IS  Prth 

PRINT  “Station: “ ;Ad( 12 >; ”  GMT  -  " ; Ad< 1  > ;  “  Date: * ;Moi "/’ j Day Yr 
PRINT  “Latitude: "iAd< 2 Lonqi tude : " ; Ad< 3 > ♦ 100 
PRINTCR  IS  1 
Rtp*«“N" 

GOTO  1758 

DISr  "PLOTTER  NOT  ON  LINC  Ml  " 

UAIT  1 
OEET  500,1 
DEEP  250.2 
GCLEAR 

OFF  TIMEOUT  llolb 

ASSIGN  0Fat  TO  •  '  close  file 

ASSIGN  80lr  TO  •  <  close  directory 

orr  KEY  4 

OFF  KEY  3 

RCTURN  !  to  menu 


0R0UN0(Depth,4 ) 


1  close  file 
!  close  directory 


!  to  menu 


tartd:  ! 1  *  Pr soar as  f or  Disk  i npu t  *  *  *  * 

. . . 

INPUT  Insert  data  disk  and  press  continue",'/* 
ASSIGN  80 ir  TO  "0IR"&Msu* 

ENTER  80ir;0ir<»)  i  read  directory  file 
PRINT  Pa* 

PRINT  “Last  AXBT  file  -  ";0ir<3)i"  out  of  8" 
PRINT  "Last  AXBT  entry  »  ";Dir(4)s"  out  of  10" 
WAIT  2 


V  V  *  *\. «  A  #  • 


,•  %  %  /• 
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1775  '  IP  Dirt  4  )-0  THEN  1777 

1776  60T0  1779 

1777  Oir<4)-0  !  reset  axbt  entry 

1778  Olr( 1 >-0ir( 1  )+l  I  uses  next  file 

1779  Dir<3)-0ir< 1 )  !  lest  exbt  file  used 

1780  0ir(4  )-0ir(4 >+l  !  update  exbt  entry 

1781  IF  0ir( 4  )> 18  THEN  1777 

1782  IF  Olr <  3  )<9  THEN  1787  !  last  exbt  file 

1783  BEEP  I  if  disk,  is  full 

1784  PRINT  “THIS  OISK  IS  FULL  *««*••••* . . 

1 /bb  A3oit3i.'  iSUii  TO  » 

1786  GOTO  Startd 

1787  Fn-0ir(3)  !  file  number 

1788  Jb«Fn«20 

1789  ASSIGN  ©Fat  TO  *FAT“&UALS< Fn  )&MsuS  I  open  file 

1790  RETURN 

1791  . . . . •••••••• 

1792  Realset:  !  prepares  for  real  time  plot  *•* 

1793  . . * . . 

1794  RtpS-“T“  I  flao  for  real  time  plot 

1795  GOSUB  Plot  •  oo  to  olottinQ  routine 

1796  PENUP 

1797  MOVE  Tm in. Dm in 

1798  RETURN 

1799  . . . . 

1800  At--  !•••*•  freauency  to  temperature  conversion  ..... 

1801  . . 

1802  Depth-Bx*T ime  !  depth  eauation 

1803  IF  Data-0  THEN  Oata-1 

1804  Temp-9000000/Data 

1805  Tempi *-126. 662+. 2 1995* Temp  -. 0001 70596 .Temp" 2+7. 70534E-8.Temp*3 

1806  Temp*Templ-1.79581E-ll.Temp''4U.73823E-15.TempA5 

1807  RETURN 

1808  . . . . . 

1809  Save:  !••  stores  data  in  das  format  ............... 

1810  . . . 

1811  PRINT  "Storina  data" 

1812  Ent-Oirt  4  )-i  ! CURRENT  XBF  ENTERY  -  1 

1813  Dirt Jb  )-Ad< 18)/100  (DATE  YY.MMDD 

1814  Dirt Jb+1 >-Adt 1 >  (TIME  HHMM.SS 

1815  Dirt  Jb+5  >« 1  !  i-DATA  TYPE  FOR  XBT 

1816  Dirt Jb+8 >«0ir(0p  +  i  >  !  LAST  XBT  IN  FILE 

1817  FOR  1-1  TO  Nots 

1818  IF  Ts<I><8  THEN  1820 

1819  NEXT  I 

1820  Adt  4  )-0pt ( 1  )  !  depth  of  8C 

1821  Adtll)-Ts<5)  !  surface  temp 

1822  GOSUB  Shrt 

1823  OUTPUT  ©Fat. Ent *20+ 1 ;0ut t •  > 

1824  FOR  1-2  TO  20 

1825  FOR  J-l  TO  36 

1826  IadtJ>-INT(Tst<I-2)»36+J >• 100  > 

1827  NEXT  J 

1828  OUTPUT  ©Fat, Ent»20+I; Iadt » > 


1829  NEXT  I 

1830  ASSIGN  0Oir  TO  ‘DIR'&Msu* 

1831  OUTPUT  Mir;Dir(  •  ) 

1832  PRINT  L2S,  “File  number* ;0ir< 3 > 

1833  PRINT  ‘Entry  number  *;0ir<4),Llt 

1834  BEEP  200.1 

1835  RETURN 

1861. . . . . . . . 

1862  Shinit:  . . . . . 

1863  . . . . * . . 

1864  !  CREATE  LOOKUP  TABLES  FOR  SHORT  BCD  FORMAT 

1865  ! 

1866  DIM  Cv$(20].Conb(2S6  > 

1867  INTEGER  Conl( 100  ).Con2( 100  ) 

1868  IF  Intxx«l  THEN  1711 

1869  FOR  Jt-0  TO  99 

1870  OUTPUT  Cv*  USING 

1871  Coni ( J I  +  1 >«SHIFT<  SHIFT ( UAL <  CvSl 1 . 1 T  > , -4  )+VAL<  Cv$l 2 .2  J ) . -8  ) 

1872  Con2<  J.lH  )”SHIFT (  UAL(  Cv$[  1 , 1 J  > ,  -4  M-UAL(  CvSt  2,21) 

1873  NEXT  J.1 

1874  FOR  J.i-1  TO  100 

1875  Conb<Con2< J i >+l >-J i-1 

1876  NEXT  J1 

1877  RETURN 

1878  . . . 

1879  Shrt:  . . . . 

1880  . . . . ****•• . . . 

1801  i  EMULATE  HP984S  SHORT  DCO  FORMAT  AND 

1882  !  OUTPUT  INTO  INTEGER  ARRAY  Out 

1883  ! 

1884  TOR  J i-1  TO  18 

1085  OUTTUT  Cv*  USING  “# ,  S.OOODDDE" ;  AdC  J.i  ) 

103G  Siaq=0 

1887  Esia-0 

1080  IT  Cv*UO.101-‘  ”  THEN  Es ia=  -32768 

1 3C3  IT  Cv*l 1,11*”-"  THEN  Giaa-256 

1330  IT  EsiaCO  THEN  1093 

1031  Out  (  J.i .  I  )“SHIFT<  <  UALC  Cv*[  11.121)1  >.  9  UE8ia»SiQa*Con2<  HUAL(Cv*I3.4  J ) 

1332  GOTO  1034 

1333  Out( J i .1  )=GHirT ( G3  < VALC  Cvtl 11,12  3)),  9 > tEeiatSiaa tCon2( ltUAL(Cv*C3.4] 

1094  Out(  J 4.2 )»Conl(UAL(Cv*[5.G)  )»1 >+Con2< VALC Cv*[ 7,8 1  >♦ 1 > 

1895  NEXT  Ji 

1G3G  RETURN 

1337  . . . . *#• . 

1898  Ushrt:  |ui»Mmni.niini*niiuiuiim 

1033  . . . . . . .  ••••«•«• 

1900  i  DECODE  HP9045  SHORT  CCO  NUMCCR5  READ 

1301  !  REAO  r ROM  INTCGCR  ARRAY  OutC 10.21  AND 

1302  '  LOAD  INTO  Adt 10  > 

1303  i 

1904  TOR  Ji-l  TO  10 

1305  Easoh^BI T ( Out ( J 1 . 1 > , 15 ) 

130G  IF  Exsan-l  THEN  EaOP-10  < -G3»BINANO<  SMITTIOutt  J.t ,  1  ).9  >,63  > ) 

1307  IT  Essan-O  THCN  EaOO-10 “ ( QINANDC  GIUTTC Out <  J.i .  I  > . 3  > .63  > H  > 

1900  Ad<  J  1  )-Conb<  SHIP T<  Out  <  J  i  ,2  l.OHl  > •  .000H Conb< DINANDC Out (  J.i ,2  >,255  > *  1 

1303  AdC  J  i  >«<  1  2  >8 1 TC  Out  (  J  f .  1  )  .0  >  )'E>.pp»  <  Conb(OINAND(Out(  J.i ,  1  >,2SS  >+l  >*.0 

1910  NEXT  J| 

1311  RETURN 
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